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Abstract 

The generating function for pjv(n), the number of partitions of n into at most N parts, may be written 
as a product of N factors. We find the behavior of coefficients in the partial fraction decomposition of this 
product as N — > oo by applying the saddle-point method, where the saddle-point we need is associated to 
a zero of the analytically continued dilogarithm. Our main result disproves a conjecture of Rademacher. 
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1 Introduction 

1.1 Rademacher's coefficients 

Let p(n) denote the number of partitions of n. The generating function for p(n) is an infinite product and 
Rademacher, in HRad73l pp. 292 - 302], obtained a decomposition for it 


Y /P (n)q n = l[ T —= J 2 E 


n —0 


j=i 


0 ^.h<k £■ 
(h,k )=1 


Chu{ oo) 

/ -j _ g27r ih/k^l 


(Id < i) 


(i.i) 


by using his famous exact formula for pin). The coefficients Chktipo) are given explicitly in !Rad73l Eq. 
(130.6)] with, for example. 


8 12 v / 3 


Coi 2 (oo) — -. + 


144 5616x/3 


1225 42875tt 


In this notation. Con is the coefficient of (q— l) -1 and C 012 the coefficient of ( q — 1) 2 . Truncating the infinite 
product in dl.H gives the generating function for pjv(n), the number of partitions of n into at most N parts, 
and its partial fraction decomposition may be written as 


N 


fiV/feJ 


J2pN(n)q n = ]J 


- E E 


Chu{N) 


n =0 


1 -qi ^ ^ [q-e 2 ™h/k\C 

j =1 y 0 ^h<k^N 1=1 ' 

(h,fe)=l 


Comparing (11.11 and dl .2b . Rademacher conjectured in !Rad73l p. 302] that 

Chki(N) Chktipo) as N 


00 . 


( 1 . 2 ) 


(1.3) 


Investigations in IAnd03l . IDG02L l Mun08l were inconclusive, but Sills and Zeilberger in I ISZ13I developed 
recursive formulas for Chki{N) and gave convincing numerical evidence that Chkt{N) Chkt{ 00 ). They 
saw that the points (A, Cqu{N)) start to trace curves oscillating with periods approaching 32 and with 
amplitude growing exponentially. Their conjecture I SZ13. Conj. 2.1] is that this is the true description. 

In 1 Q'S151 we found relatively simple, explicit formulas for Rademacher's coefficients Chkt{N), linking 
them to formulas of Sylvester [Syl82) and Glaisher I Gla09l . For example IO'S151 Eq. (2.12)] is 


Coit(N) = 


(-1 r(f-i)! 

IV! 


E 

jo+jl+ja-\ -h jN=N—t 


l + jo ] B jl B j 2 ---B jN l jl 2 32 ■ • • N jN 
£ J (£ — 1 + jo)! jiW ■ ■ ■ jjv! 


where Bj is the jth Bernoulli number and { ” } the Stirling number denoting the number of ways to partition 
a set of size n into m non-empty subsets. Also in lO'S15l , based on an earlier stage of the work in this paper, 
the exact asymptotic behavior of Con (IV) was conjectured. This requires the solution wq ~ 0.916 — 0.182? to 


Li 2 (w;) — 27rilog(w;) = 0 


(1.4) 


where I.i 2 denotes the dilogarithm. (It may be seen that wq is a zero of the dilogarithm on a non-principal 
branch, see Section (231 ) Set zq := 1 + log(l — wo)/(2m) so that 


w 0 = 1 - e 2 ™* 0 , 1/2 < Re(z 0 ) < 3/2. 

Conjecture 1.1. il Q'S15\ Sect. 61 We hav$\ 


Con (N) = Re 


(-2z 0 e-^°)- 


,-N 


N 2 


O 


bo | 


-JV 


N 3 


(1.5) 


( 1 . 6 ) 


^his statement is equivalent to Conjecture 6.2 in lQ'S15l where zo and wo are replaced by their conjugates. 
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(1.7) 


Equivalently, we may present (11.6b more explicitly as 


Con (IV) 


^(asin(/3 + llV)+o(l)) 


for 


U :=- log |w 0 1 ~ 0.0680762, V := arg(l/w 0 ) « 0.196576 (1.8) 


and a := | — 2izoe~ mz °\ ~ 5.39532, /3 := arg(— 2izoe~ 7rlz °) ~ 1.21367. This implies the period of Con(lV) is 
2n/V ~ 31.9631. As we will see, the numbers W(, and zq control the asymptotics for all of the Rademacher 
coefficients that we examine. 


1.2 Main results 

Write the Farey fractions of order N in [0,1) as 

:= j/i/fc : 1 ^ k ^ N, 0 ^ h < k, ( h,k ) = l|. (1.9) 

Our first result is a kind of averaged version of Conjecture !! .1 i with Con (N) replaced by 

C 011 (N) + C 121 (N) + ■ ■ ■ + C {99)(100)1 (N) = J2 C hki(N). 

h/kG^ioo 


Theorem 1.2. For an absolute implied constant 


J2 C hkl {N) = Re 
h/ fcG«^ioo 


(-2z 0 e-™°) 


w~ N 

— 7TIZQ \ ^0 


N 2 


+ o 




( 1 . 10 ) 


This has the following consequence. 

Corollary 1.3. There exists a pair {h,k) with h < k < 100 such that limjv-^oo Chki(N) does not exist. Hence 
Rademacher's conjecture that Chki(N) Ch k e(oo) as N —> oo is false. 

Proof. Expressing the right side of fll.lOt as in G3, we see that this side cannot have a limit as N —> oo since 
a f 0, U >0 and B + VN comes within 1/10 of n/2, say, infinitely often since V <1/5. But the left side 
of dl.lOb is a finite sum, so Rademacher's conjecture implies that its limit as N —> oo exists. The corollary 
follows. □ 

Theorem 1 1.2 1 is the i = m = 1 case of the next result where we extend the right side of (II .10b to include 
the first m terms of the asymptotic expansion and generalize Cqi i (A - ) to Cqu(N). 

Theorem 1.4. There are explicit coefficients c^o, c^i,... so that 


t 

Cou(N) + J2 Y,^ ih/k 

0<h/k£j?ioo j —1 


i y- j c hkj (N) 


= Re 


w. 


-N 


N e+1 


(cifi 


Cl, 1 

w 


H-h 


Q,m-1 
jym—l 


+ o 


( Kr” \ 

l ]\J-e+m +1 J 


( 1 . 11 ) 


where c^o = — 2z$e 7rzz °(27r*zo) f 1 an d the implied constant depends only on £ and m. 

See ( 15.41b for the next coefficient cyi. The reason we need to include the sum over 0 < h/k £ l^ioo on 
the left of dl.llb is given in Remark 1 6.3 1 Numerically, this sum looks to be much smaller than Cqu(N), so it 
is natural to generalize Conjecture ll.ll to: 

Conjecture 1.5. For the coefficients a t o, Q,i,... ofTheorem \1.4\ 


C 0U (N) = Re 


0*1 ( r 
N e +1 V f ’ c 


ce,i 

W 


C£,m— l\ 

+ TV" 1 ” 1 / 


+ 0 


f kol-^ A 

jyt+m+l J ■ 


( 1 . 12 ) 
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Numerical evidence for Coniecture ll.5l is given in Section[ 6 ]and the asymptotics of the next cases, C 121 (IV) 
and Ci 3 i( 7 V), are also discussed there. The following subsection outlines how the proofs of Theorems 11.21 
and ll.4l are constructed. 

As this work was being completed, I was contacted by Drmota and Gerhold who provided me with their 
paper HDG14I . They have given an independent disproof of Rademacher's conjecture by combining a Mellin 
transform approach with the saddle-point method to obtain the asymptotics of Cou{N). Their main result 
is equivalent to Conjecture II .51 in the case m = 1 though with a weaker error term. Combining their tech¬ 
niques with ours should lead to improved asymptotics and a better understanding of all the Rademacher 
coefficients. 

1.3 Method of proof 

We have from lO'S15l Eq. (2.1)] that 

p 2ixiz ( Jl-niz p 2nih/k\£—l 

ChkAN) = 2m Res -- , .. - 5 -^r- / -(1.13) 

v ’ z=h/k( 1 - e 2mlz )(l - e 2m2z ) ■ ■ ■ (1 ~ e 2niNz ) 

The right of tl. 13b may be expressed in terms of the simpler function 

^2ivi(TZ 

(j) . ^ _ g2ni2z^ . . . ^ z'j (1*14) 

and we write 

Qhka(N) := 2iri Res Q(z-,N,a). (1.15) 

z=h/k 

Expanding the numerator on the right of dl.l3b then produces 

ChkAN) = E P M (-e 2 ” ih / k y-°Q hka (N). (1.16) 

<7=1 ' ' 

The numbers Qhka{N) are slightly easier to work with than ChkAN), though of course for l = 1 we have 
ChkAN) = Qhki(N). Each QhkAN) is a component of the Sylvester wave W k , as described in Section [6721 
and expressions such as Q(z; N, a) and QhkAN) appear when counting lattice points in a polytope dilated 
by a factor — a > 0, see [ BDR02 . Thm. 1] and IBGK01I . We may also invert (II .16b to get 

QhkAN) = E ( a ~l)(e 2 ^ k r~ e C hk AN). (1.17) 

t—i k 7 

As a function of z, Q(z\ N. a) is meromorphic of period 1 when er £ Z. Fixing a positive integer a and 
summing all the residues then leads to the key identity on which Theorem ll.4l is based: 

E QhkAN) = 0 for N(N + 1)/2 > a. (1.18) 

h/k^^N 

There is a large contribution to the left of (11.181 from QoiAN) as well as other QhkAN) with k small, 
corresponding to high-order poles of Q(z; N, a). Balancing that are contributions from coefficients Qhka{N) 
with k large, corresponding to simple poles. Put 

A(N) := { h/k : y < k < N, h = 1 or h = k - l} C & N (1.19) 

and decompose 111 . 18b into 

E QhkAN) + E QhkAN) + E QhkAN) = o. 

h/k£ c^ioo h/k^^N — (^’iooU.A(iV)) h/k£A(N) 

The reason we focus on the subset A(N) is given in the next section, but it may already be noticed that, 
numerically, 

CqiAN) ~-AAN,1) as N oo (1.20) 
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for 


( 1 . 21 ) 


Ai(N,a):= J2 QhkAN). 

h/keA(N) 

Computing the residues of the simple poles lets us describe lll.21b more explicitly as 


„ , T v- 2(—l) fc fm 

A\ ( N , a) = Im 2_^ —P5— ex P ( ' 


k 2 


-<k^N 


\ 2 


-N 2 -N + 4cr 


3N 


rtiv—fe (i/^) 


( 1 . 22 ) 


for a £ Z, where we write 

m 

Il m W ~l[2sm(TrjO) (1.23) 

i=i 

with IIoW := 1/ following Sudler's notation in [ Sud64 | except that we don't take the absolute value. The 
main part of the proof of Theorem [L4] then consists of establishing the following two results. Recall w {} and 
zq from US. 

Theorem 1.6. With bo = 2zoe~ 7rzz ° and explicit b±(a), b 2 {cr), ■ ■ ■ depending on cr £ Z we have 


Ai(N, a) = Re 


—N 


N 2 



jv™-i J + \ N m + 2 J 


for an implied constant depending only on a and m. 

Theorem 1.7. There exists W <U := — log |u>o| ~ 0.068076 so that 

J2 Qhk*(N) = O (e WN ) 

h/ke^N — (^iooOA(N)) 

for an implied constant depending only on a. We may take W = 0.055. 

The proof of Theorem 11.61 is carried out as follows. In Section [2] we derive the sum for A-\ (N. a) in 
dl.22b and also give some results on the dilogarithm we will need. Section[3]is quite technical and includes 
estimates of the sine product IlAr-fcC^/^) using Euler-Maclaurin summation, where the number of terms 
required is proportional to k and N. The sum Ai ( N, a) is replaced by an integral in Section[4] and in Section 
[5] the saddle-point method is introduced and applied. The required saddle-point is zq and with work of 
Wojdylo |Woj06| we explicitly get the full asymptotic expansion of A\ (N. a). This proves Theorem |1.61 

See Section[ 6 ]for a summary of the proof of Theorem ll.7l The bounds required for Qhka(N) in this proof 
can be obtained directly in most cases, but three families also require saddle-point arguments, with these 
saddle-points corresponding to further zeros of the dilogarithm on other branches. The details are carried 
out in the companion paper IQ'Sall . 

Linear combinations of Theorems ll.6l and ll.7l then give Theorem 11.41 in Section 15.41 In Section [ 6 ] we also 
discuss extensions and generalizations of our results and applications to the restricted partition function 
and Sylvester waves. 


2 Preliminary material 

2.1 The residues of Q{z ; N, a) 

For Q(z',N, a) defined in 111.14b with a £ C, we clearly have 

Q(z-,N,<t) = Q(-z;N,<t), (2.1) 

Q(-z; N , a) = (-1 ) n Q(z; N, N{N + l)/2 - a) (2.2) 

and, if er £ Z, 

Q(z+ 1;N, a) = Q(z; N, a). (2.3) 

As a function of z, Q(z ; N. a) is meromorphic with all poles contained in Q. More precisely, the set of 
poles of Q(z; N, a) in [0,1) equals .^Ly, the Farey fractions of order N in [0,1). 
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Theorem 2.1. For N e and <r e Z 


-pAr(-cr) 


27ri ^ Res Q(^; AT, a) = < 0 


h/k^^N 

Proof. We have 


z=h/k 


if cr ^ 0 

if 0 < cr < iV(iV + l)/2 (2.4) 


(-l) Jv p JV ((7-iV(iV+ 1 )/2) zf A(A + l)/2 < a. 


^PAr(n)( 


2-Kinz 


n =0 


N 

n 

7 = 1 


1 


1 — e 2 ,r h 2 


(2.5) 


and since pn{ji) ^ p(n ) sC 2 ra_1 , the number of ordered partitions of n, we see the left side of (12.51) is 
absolutely convergent for Im(z) large enough. (Better bounds for pn{u), p(n) imply absolute convergence 
for Im(z) > 0. See for example IIPri09l . employing the dilogarithm.) Hence, for lin(u') large enough. 


W -\-1 


Q(z; N, a) dz = 


| 0 if a > 0 

\p N (-cr) if t7<0. 


( 2 . 6 ) 


Let 1Z be the rectangle in C with upper corners w, w + 1 and lower corners v, v + 1 where liu(i') < 0. 
Integrating around R in a positive direction and choosing Re(w) = Re(u) between 0 and the next pole to 
the left, 

~ (2.7) 


/ Q(z; N,a) dz = 2iri Res Q(z;N,a). 

Jn , z=h/k 


h/k^i^N 

The integral along the top of R is — 1 times j 2.6 b . The integral along the bottom can be made arbitrarily small 
by letting Im(i>) —> — oo provided cr < N(N +1)/2 and the integrals along the vertical sides cancel with 12.3b . 
If cr ^ N(N + l)/2 then use 12.2b . This completes the proof. □ 

Theorem 12.11 for negative integer a is a restatement of a special case of Sylvester's Theorem. See for 
example I IO'S151 Sect. 4]. 

Each h/k € RFn is a pole of Q(z; N, a) of order s = \_N/k\ . Equivalently, h/k is a pole of order s exactly 
when 

<k^~. (2.8) 

Thus 2niQ(z; N, a) has one pole of order N in [0,1) at h/k = 0/1 with residue Qoi<r{N). The next highest 
order pole has order [TV/2J at h/k = 1/2 with residue Qi 2 a(N)- By 12.8b . h/k is a simple pole when N/2 < 
k < N and the residues of the simple poles of Q(z: N, a) may be computed quite easily. 

Proposition 2.2. For N/2 < k ^ N 


Qhka(N) = 


(-D 


fe+1 


k 2 


■ exp 


—7r ih(N 2 + N — 4a) 


2k 


x exp ( f(2Nh + N + h + k-hk )) n 2sintoW 


JV-fc 


Proof. With 11.13b . write 


QhkAN) = Res 


27t ie 


2n icrz 


z=h/k [(1 - e 2viz ) •••(! — e 2 -rri(fc—!)«)] (1 _ e 27rifcz) ^ _ £ 2ni(k+l)z ) ...(!_ e 2niNz )j 


Then 


Also 


1 


Res , 

z—h/k 1 — 6 


2nikz 


-l 

2nik 


(1 — 0(1 — C 2 ) •••(! — C fe_1 ) = 


(2.9) 
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for { = e 2 " th / k ‘ a primitive fcth root of unity, by 60'S151 . Lemma 4.4] for example. Hence 


Qhk„(N) 


""fc 5 


N 


n 

j=k +1 


1 

^ _ glirijh/k 


^l-niah/k 

""fc 5 


N—k 


n 


i 

^ _ g2itijh/k ' 


A straightforward calculation, with 0 ^ m < k, shows 

TT 1 h. 1 

1JL 1 - = “» {— P - t (m + 1( j J h 2 sin Mh/k) 

and combining this with (12.10b and simplifying completes the proof. 


( 2 . 10 ) 


( 2 . 11 ) 

□ 


2.2 Products of sines 

Recall our notation 11 .231 . Then for integers k > h ^ I with (h, k) = 1 


Yl m \h/k) = n 


2 sin(7rj/i/fc) 


(0 < to < k) 


is a real number. For example, 

Y[- k - 1 {h/k) = {-li h - 1 ^' 2 \ 

follows from setting m = k — 1 in 12.lit and using 12.9t . 

With Proposition 12.21 we see that 


\Qhka(N)\ 


n 


-i 

N-k 


(. h/k) /k 2 


(A/2 < k < N, a G R) 


( 2 . 12 ) 


(2.13) 


and so the size of Qhka (A) is controlled by the sine product [jj (h/k) for to = A — k. As to varies we need 
to know how big JX m (h/k) can be. For example, if h = 1 and k is large then the first terms 

1 k 1 k 

2sin(7rl/fc) 27 t’ 2sin(7r2/fc) 47r’ 

are all greater than 1. The maximum is reached with 


1 


2sin(7rl/fc) 2sm(n2/k) 


2sm(n(k/6)/k) k ^ 


rifc/6 


(2.14) 


since after that the factors become less than 1. If h = 2, the maximum value of Jim (2/fc) is reached for 
to = k /12 and this value is approximately the square root of 12.14b . For other values of h the maximum of 
the product does not become as large because values greater than 1 are multiplied by more values less than 
1. An exception is when h = (k — l)/2 since here again large products can build up. We illustrate this with 
Figure [l] which graphs 


'Ll h/k) := max 

0 ^m<k 


iog|n m (v*o| 


(2.15) 


for 1 ^ h < k with k prime and equalling 211. We see that the largest values of 4 /(h/k) are for h € 
{1,2, (fc — l)/2} (and symmetrically h £ {k — l,k — 2,(/c + l)/2}) with exactly these values greater than 
U ~ 0.068076. These observations are made precise in Section [All 

So, among the simple poles of Q(z; A, a), Figure [H leads us to expect that the largest contribution to the 
left of (11.18b should be from the sum .Ai(A, a) as defined in (11.21b . 

With (12. lb and (12. 3b we obtain the identity 


27t* Res Q(z; A, a) = 2iri Res Q(z;N,a ), (a£Z). 

z=l — h/k z—h / k 
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Therefore, assuming a G Z from now on. 


Ai(N, a) := 2ni Res Q(z\N,a) 

^ J z=h/k 


h/keA(N) 


= 2 Re 


27t* Res Q(z;N,a ) 

z=l/k 

%-<k<N 


So setting h = 1 in Proposition 12.21 a nd simplifying yields 111. 22b . 


2.3 The dilogarithm 


We assemble here some of the properties of the dilogarithm we will need. See for example IMax03 l. |Zag07| 
for more details. Initially defined as 


OO yi 

Li 2 (z) := V ^2 for M < (2-16) 

rt=l 

the dilogarithm has an analytic continuation given by 

- [ log(l -u)— (2.17) 

JC(z) u 

where the contour of integration C(z) is a path from 0 to z £ C. This makes the dilogarithm a multi-valued 
holomorphic function with a branch points at 1, oo (and off the principal branch another branch point at 0). 
We let Li 2 (z) denote the dilogarithm on its principal branch so that Li 2 (z) is a single-valued holomorphic 
function on C — [1, oo). 

To see why the dilogarithm appears in our calculations, recall that Qikcr(N) = — k~ 2 q a TT/li 1/(1 — <I :i ) 
for q = e 27rj / fe and to = N — k by ( 12.10b when N/2 < k ^ N. Then for 0 ^ to < k 

n = ex P (y E lo s i 1 - Q j )J 

" eXP (" Io l0g (1 “ qX) dX ) = 6XP (“2^i / l0g (1 “ z) t) 

= exp ^A(Li 2 ( e 2--A) _ Li 2 (l))^ . (2.18) 

Of course, the approximation "cs" above is not very accurate. We make it precise by separating the argument 
of TP" 1 1/(1 — o- 7 ) from its absolute value TT 1 (1 /fc ) . as in <12.1 ltl. and then carefullv estimatine IT/ \l/k) 
in Propositions 3.9 and 3.14. The dilogarithms in (12.181 1 will reappear in <15. lb and (15.2b . 
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For z £ C we have the functional equations 


Li 2 (l/z) = -Li 2 (z) - Li 2 (l) - ^log 2 (-z) 


zt\ 0,oo), (2.19) 

Li 2 (l - z) = - Li 2 (z) + Li 2 (l) - log(z) log(l - z) z g (- 00 , 0] U [1, oo) (2.20) 


from IMax03i Eqs. (3.2), (3.3)], where we mean the principal branch of the logarithm on C — (—oo, 0]. Re¬ 
placing 2 by e 27riz in (12.19b and 12.20ll gives: 


• For m £ Z and m < Re(z) < m + 1 

U 2 (e~ 27riz ) = - Li 2 (e 27riz ) + 2 tt 2 (z 2 - ( 2 m + 1)2 + m 2 + m + 1 / 6 ) . 


( 2 . 21 ) 


Let (—too, m] denote the vertical line in C made up of all points with real part m £ Z and imaginary 
part at most 0. Then for m— 1/2 < Re( 2 ) < m + 1/2 and 2 qL (—i 00 , m] 


Li 2 (e 2viz ) = - U 2 (1 - e 2 ™ z ) + Li 2 (l) - 2 t n(z - m) log (l - e 2viz ). 


We may describe Li 2 ( 2 ) for 2 on the unit circle as 


^ cos(27ma:) ,, , ... 

Re(Li 2 (e 2 )) = - 2 - = n 2 B 2 (x — [x\) (x £ 


71=1 


t /t . / 2ir?r\\ sin(27r nx) N 

Im(Li 2 (e 2 )) = 2^ —-- = Cl 2 (27ra;) 


(x £ 


71=1 


where B 2 (x) := x 2 - x + 1 /6 is the second Bernoulli polynomial and 

C1 2 (0) := — f log |2sin(a;/2)| dx (6 £ 
J 0 

is Clausen's integral. Note that Li 2 (l) = £(2) = 7r 2 /6. 



( 2 . 22 ) 

(2.23) 

(2.24) 


(2.25) 


The graph of C1 2 (0) is shown in Figure |2] It is an odd function with period 27t and maximum value 
C1 2 (7t/3) ~ 1.0149416 indicated. Since, with 112.2411 . 

C1 2 (0) = (Li 2 (e ie )-Li 2 (e 7e )) /(2 i) 

= {U 2 (e l9 )-U 2 (e-' e )) /(2i), (2.26) 

we may use (12.261 1 to obtain the analytic continuation of C1 2 (0) to 6 £ C with 0 < Re (ft) < 2-k for example. 
Another approach to this continuation combines (12.231 1 and 12.2411 to get, for m^z^m+l 7 m£Z, 

Cl 2 (27 tz) = — iLi 2 (e 27 ” z ) + itt 2 (z 2 — (2m + 1 )z + m 2 + m + 1/6) . (2.27) 

Then the right of 12.2711 gives the continuation of Cl 2 (27tz) to z £ C with m < Re(z) < m + 1. 

As z crosses the branch cuts the dilogarithm enters new branches. From I Max03 l Sect. 3], the value of 
the analytically continued dilogarithm is always given by 

Li 2 (z) + 4 tt 2 A + 2niB log (z) (2.28) 
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for some A, B £ Z. For example, the simplest way to get to a branch corresponding to J2.28b is to circle 
z = 1 once in the negative direction, then circle z = 0 in the positive direction A times and circle z = 1 in the 
negative direction B — 1 times (using the opposite directions it A or B is negative). 

Zeros of the analytically continued dilogarithm will play a key role in our asymptotic calculations and 
in lO'Sbfl we have made a study of them all. When the continued dilogarithm takes the form ll2.28b with 
B = 0, there will be a zero if and only if A + 0 and, for each such A, the zero will be unique and lie on the 
real line. The cases we will require have B ^ 0. In these cases there are no real zeros so we may avoid the 
branch cuts and look for solutions to 


Li 2 (z) + 47 t 2 A + 27 riB log (z) = 0 (z£C, 2 ^ (— oo, 0] U [1, oo), A,B £ Z). (2.29) 

The next result is shown in Theorems 1.1 and 1.3 of lO'Sbl . 

Theorem 2.3. For nonzero B £ Z, (12.29b has solutions if and only if—\B\/2 < A < \B\/2. For such a pair A , B the 
solution z is unique. This unique solution, w(A, B), may be found to arbitrary precision using Newton's method. 

Sketch of proof. By considering the two curves where the real part and the imaginary part of (12.28b vanish, 
it can be shown that they intersect if and only if — \B\/2 < A A \B\/2. It can also be shown that w( A. B ) is 
close to e 2nlA / B and from this starting point (the case A = 0 needs an adjustment) Newton's method will 
always converge to w(A, B). □ 

By conjugating (12.29b it is clear that 

w(A, —B ) = w(A, B). 


So for nonzero B the first zeros are w(0, 1) and its conjugate w(0, —1). We have 

w(0, —1) w 0.9161978162 - 0.18245889721 
and this zero was denoted by w o in Section [l] The next few zeros are 


u>(0, -2) « 0.9684820460 - 0.10953110651 

w(l, -2) « -0.9943069304- 0.06488893181 
w(- 1, -3) w -0.5459030969+0.88123074231 
w(0, -3) Ri 0.9832603795 - 0.07775963891 
to(l,-3) « -0.4594734813- 0.84853503801 


where w(0, —2) and tv (1, —3) will be required in Sectional 
Define 

_Li 2 (e 2 -)+Li 2 (l)+47r 2 d 

Pd(z) ~ ----, (2-30) 

Imz 

a single-valued holomorphic function away from the branch cuts (—loo, n] for n £ Z. In Section[5]we will 
require the solution of p' 0 {z) = 0 (and in IIP Sail solutions to p' d (z) = 0 more generally). 


Theorem 2.4. Fix integers m and d with —|to|/2 < d < |m|/2. Then there is a unique solution to p' d (z) = 0 for 
z £ C with m — 1/2 < Re(z) < to + 1/2 and z yt (—loo, to]. Denoting this solution by z* , it is given by 


and satisfies 


* i 

Z = TO + 


log(l — w(d, —TO.)) 


27t1 


Pd{z*) = log (w(d,-m)). 


Proof. Note that 

^ Li 2 (e 2mz ) = -2t ri log (l - e 2niz ) 
for z not on any of the vertical lines (—loo, n], n £ Z. So 


(2.31) 

(2.32) 


Pd{z) 


zp'd(z) = -r {zPd(z)) 
az 


Li 2 (l) + 47r 2 d Li 2 (e 2 


27ri 


27t1 


= log (1 - e 2 ^) 
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and hence 


p' d (z) = ~ (p d (z) - log (1 - e 27viz )) . (2.33) 

Similarly 

1 / 27T? • P 2viz \ 

= -- (VaW + Y^r) ■ (2.34) 

With d2.30t we may expand (12.33b into 

2mz 2 p' d (z) = Li 2 (e 2mz ) — Li 2 (l) — 47r 2 d + 27U;zlog (l — e 2mz ) . (2.35) 

Applying the functional equation d2.22b to d2.35b implies 

27 jiz 2 p' d (z) = — Li 2 (l — e 2mz ) — 4TT 2 d + 2mm log (l — e 2mz ) 

for m — 1/2 < \\q(z) < m + 1/2. Letting w = 1 — e 2niz , we are now looking for solutions to the equation 

Li 2 (w) + 47r 2 d — 2mm log (w) = 0 (2.36) 

and Theorem 12.31 gives the unique solution as w(d. —m) when — \m\/2 < d sC \m\/2. The formula d2.31b 
follows and then d2.33b implies d2.32b . □ 


3 Estimates for the sine product Yl^ih/k) 

Sudler in | [Sud64j approximates J~[ m (0) us i n g the first term of the Euler-Maclaurin summation formula and 
finds that the 9 £ (0,1) that maximizes | ]/[ TO (0)| i s approximately xo/m where xq ~ 0.791227 is the x value in 
(0, 1) where /t. (C1 2 ( 2 tt x )/( 2 ?r x)) vanishes. Wright in IWri641 uses more terms in the summation to get more 
detailed results, as do Freiman and Halberstam in [ FH881 . We use similar techniques in the next subsection 
but require arbitrarily many terms of the summation formula. 

3.1 Euler-Maclaurin summation 

We need to estimate the size of ]~[ (0) accurately and also replace it with a continuous (and later holomor- 
phic) function of to. 

Let p(z) := log ((sin z)/z), a holomorphic function for \z\ < n that satisfies p(—z) = p(z). Also 

cot(7T2) = — + p'(irz) (3.1) 

7 TZ 

and so 

cot (fe) M = [ n ^k+i + P ik+1) (nz)- (3-2) 

Proposition 3.1. For to, L e Z^i and —1/m <9 < 1/m with 9 f Oive have 


IU*) = (|0| 


9 \ m f 2sm(mn9)\ 1 ^ 2 ( C 1 2 ( 27 tto 0 ) 


exp 


2t t0 
/ L—l 


x ex P 1 cot* 2 * 2) (mn9)J exp(T L (m,9)) (3.3) 


for 


m /■” - |ij) , 31 „ , r B 2 L -B 2 L (x-[xl) , 

T l( ™,«):=(»«) ^ - — -/ 'H)d<+| 2 L(x + „,P dx ~ 
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Proof. Write 


n„,w 


(2rt) m m!fl 

A A 7T7 6/ 

i=i 

m m 

(27r0) m m! exp (p(nj9)) = (27r6>) m m! exp p(j:j9). 
3=1 3=1 


(3-4) 


With Euler-Maclaurin summation, as in |Rad73. Chap. 2] or HQlv74l p. 285], we obtain for \0\ < 1/m, 


"<■ pm | 

p(njO) = / p(ttx8) dx + — (p(7rm0) — p(tt09)) 
3=1 J ° 2 


+ 1 {p (2£ ^(Trm^) - p (2t 1} (tt06»)| + R L (m, 8) 


(3-5) 


where L > 1 and 


f?i(m, 8) := (ttO) 2L f 

Jo 


2L B 2 l — B 2 l(x — [a:]) 


(2 L)\ 


P^ 2L \ttx8) dx. 


(3-6) 


The integral in (13.51) may be evaluated using (12.25b to get 


and therefore 


ruw = 


9e 

\9\m 


pm 

/ p(ttx9) dx = — to log |27 tto0| + m — 

J o 


sin(7rm.0) \ / Cl2(27nn0) 


Cl2(27TTO0) 
2 tt0 


(0/0) 


7T 7710 


exp 


2nd 

(L -1 


x exp ^ ( 7rm6 ')J exp ( B L {m,9)). (3.7) 


Stirling's formula is 


L—1 


logr(m) = (m - -)logm - m + - log2^ + ^ - i^to 2 *" 1 + 5i(w) 


with 


£7, (to) := [ 

Jo 


B 2L - ^2i(a; - |xj) 


0 2 L{x + to) 


2L 


dx 


(3.8) 


(3.9) 


as in IQlv74l (4.03) p. 294]. Hence 

Gs)” m! = (s)” mr(m) = < 2 ™) 1/J ®P (E a(2<-”)m*->) “ p(Sl(m)) - < 3 - 10 > 

Putting (13.10b into ( 13.7b , recombining the two sums with ( 13.2b . and setting Tjf m, 9) := Rif m, 9) + ,57. (in) 
completes the proof. □ 


3.2 Derivatives of the cotangent 

We next examine the cotangent function and its derivatives in detail. For all z £ C, 

, . 2i 2 i 

cot(Trz) = » + e2Tfa _ - = -i - e _ 2jriz _ x ■ (3-11) 
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For k ^ 1, by induction. 


k-\-l s 

cot (fc) (z) = (—l) fc (2z) fc+1 - 1)!< 

r—l k 


k + 1 


r J (e 2iz — l) r ’ 

= (-l«-2^ +1 E(-l)!{‘t 1 }(F^^- 

As in 1GKP94 ;. Chap. 6] these Stirling numbers satisfy the relations 


k 

r—l 


k + 1 
r 


r—0 


1 T . \ x ( x - 1) • • • {x - r + 1) = 


For k ^ 0, cod *9 (tt;;) is clearly holomorphic in C except for poles when zfZ. 
Lemma 3.2. For c > 0 raid k £ Z^ 0 


J2? k e - ce - ^ k\ 


fc+i 


Proof. The result follows by comparing the series to the integral 

T(k + 1) _ ft! 


x k e cx dx = 


r k -\-1 /'•fc+l 


Theorem 3.3. For all nonzero z £ Cwith —1/2 ^ Re(,s) < 1/2 we have 


COt^(7 TZ) 

Also, for all z £ C with \y\ ^ 1, 


^ So ,fc + 20- 


k! 


fc+1 \ I ~ I fc+1 


* k+1 I e _7I ’l' y L 


cot (fe) M | < ^ e " |y| • 


Proof. Bv IRad73l 111.111 and ll3Tl) 

/! \ 2 2 ’ \B 2r \ 2r— 1 n I ^ \ 

-pW = L ( 2r )i w (l w l <7r ) 

so that all the coefficients of —p'(w) are positive. Hence, with w \ < n and the bound 

\B 2n \ - 2 




(2n)! 3(27 t) 2 " 


from »Rad73l (9.6)], we have 

OO 

\p'{w)\ < -p’(\w\) 


r—l 


Letting f(t) := 


3(1 — t/ 7r) 


we see that 




2r— 1 oo 


kl 

3 7r fc-i(l _ ^/7i-)fc+i 




7T J 3(1 — \ w \/’ k ) 


(3.12) 

(3.13) 

(3.14) 


□ 


(3.15) 


(3.16) 


(3.17) 


(3.18) 


Since the power series coefficients of / are greater than the corresponding power series coefficients of —p', 
and all coefficients are positive, it follows that the coefficients of f <k> are greater than the corresponding 
coefficients of — p^ k+1 \ Therefore 


,(fc+i) 


O) 


€ 


k! 


3n k 1 (1 — \w\/n) k+1 


(k ^ 0, |ru| < 7r). 


(3.19) 
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With d3.2b and (13.19b we have proved that 


C0t^(7Tz) 


€ 


k\ 


3 7T fc+1 \Jz | fc+1 (1 — |z |) fc+1 

Next we assume y f 0. Formulas (13.11b . (13.12b and (13.13b imply 


N < !)• 


(3.20) 


fc+i 


cot^(7Tz) < So,k + 2 fc+1 £> - l)l| k + r 1 } (1 ^ e-Uivi y 


^ So,k + 


2 k+1 


(1 — e —27r|yhfe+l 

' ' r =1 

fc+1 


fc+1 

^ r fcg-27r|y|r 


\ * 


fc+1 


^ So ’ k+ {l-e-^\y\) 


e -*\v\ y' r k e -*\y\ r 


r—1 


where we used that (r — l)l| fe ^ 1 1 r k which follows from the relation on the right of (13.14b with 
Lemma lT2l applied to (13.21b shows 


co 




^ So,k + k\ 


\ fc+i 

) {yfO). 


(3.21) 


x = r. 


(3.22) 


^7r|y|(l — e 27r l^l) / 

Then d3.15b in the statement of the theorem follows by combining d3.20b for \y\ b 0.55 and d3.22b for y ^ 
0.55. Finally, d3.16b follows from d3.22b . □ 


3.3 Initial estimates for Yl m (h/k) 

Stirling's formula implies 

2^Q"<n!<3^Q" (neZ^) 

and it follows that 

(n — 1)! < 3 {n/e) n . (3.23) 

Since |f? 2 n(£ — [^J)| ^ \B 2n \ as in Il01v74i Thm 1.1, p. 283], we see that B 2n — B 2n (x — L^J) has the same 
sign as B 2n if it's non-zero. Therefore the terms in d3.8b alternate in sign and it follows that ( TQlv741 (4.05) p. 
294]) 

\B 2L \ 


|Sz,(to)| ^ 


2L(2L - l)m 


2L-1 ‘ 


Employing d3.18b and d3.23b in d3.24b gives 


lo , 7T (2L-2)! 7T ( 2L — l 
\S L (m.)\ < — 7 “- zrn~~^r < - ' 


2L-1 


(L, m e Z^i). 


2L — 1 


6 (27t m) 2L ~ l 2 \ 2nem 

Lemma 3.4. For and —1/m < 6 < 1/m we have 

\R L {m,Q)\ < — ( 

1 K ’ 3 \27re(l — m\9\) j 

Proof. With d3.6b . d3.18b and the inequality \B 2n — B 2n (x — [xj)| < 2\B 2n \ we have 

27T 2 

w.mr , 

/ 0 


(3.24) 


(3.25) 


(3.26) 


o 2 pm 

\ R L (m,9)\ < Jo p{2L) 


dx 


3 V 2 


2L —1 


n(2i-l) 


( 7 rm$) 


Then applying d3.19b shows 


|f? L (m,0)|<-(2L-2)! 


\e\ 


27r(l — m|0|) 


2L-1 


The result follows with d3.23b . 


(3.27) 

□ 
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We now concentrate on the case where 0 = h/k for relatively prime integers k > h ^ 1. We think of h, k 
as fixed with integer m varying in the range 1 $5 rn < k/h. 

Lemma 3.5. For 1 ^ m < k/h we have 

|Ti(m,/i/fc)l < TT 2 h/ 18 +l/ 12 . (3.28) 

Proof. Note that since m is an integer 

1 ^ m < k/h =>■ h ^ mh ^ k — 1. (3.29) 

Consequently 1/(1 — mh/k) ^ k and using this in (13.27b gives a bound for R\. Bound Si with 13. 24l l. □ 

Define 

0/2 

c(h) := exp(7r 2 /i/18 + 1/6). (3.30) 

(We increased 1/12 in 13.281 to 1/6 to ensure c(h ) > 1 for h Is 1, as needed in Proposition 13.91 1 The next 
result gives us our initial estimate for JT (h/k). 

Proposition 3.6. For 1 ^ m < k/h 


n J(h/k) < c(h) exp Ch(2nmh/k)^J . 

Proof. Combining Lemma [3.51 with Proposition l3.1l shows 


Note the simple inequality 

and hence 


1 < 


€ 


1 ^ 


sin a; 2 


€ -— € 


1 / b \ 

exp(7r 2 /i/18 + 1/12) exp ( — — Ch^nmh/k) j . 

(3.31) 

7 ; ( 7t/2 ^ X < 7r/2) 

(3.32) 

7 r 


— for 0^e<a;^7r — e. 

2e 


tt/ k and so 


,s < 4 (1 < rn < k/h). 

(3.33) 


Inequalities 13.311 and 13.331 complete the proof. 

Proposition 13.61 implies that 

T& (h/k) < c(h) for k/2h < m < k/h , 


□ 


(3.34) 


since 02(0) ^ 0 for n R. 0 R. 2tt. For the rest of this subsection we focus on m in the range 1 sk ra -C k/2h. 
Our next goal is to show that for m near the end points of this range the product Wf n (h/k) is also quite 
small - see Figure |3] We first develop a simpler version of the bound in Proposition 13.61 

Lemma 3.7. For 0 < 9 < 7r we have 


C1 2 (0) < 6> — 6>log6» + 6» 3 /54, 
02(71 — 9) < 01og2. 

Proof. Integrate 13.171 twice and use 12.251 to show that, for 0 ^ 9 ^ tt, 

\B 2n \ 


ci 2 (0) = e-e\oge + Y^ 


2n(2n + 1)! 


q2ti+1 


^e-eioge + e^Y^ 


e 2 


3 ^ 2n(2n+ 1)(2 tt) 2 ”' 


(3.35) 

(3.36) 
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1 


The series above is bounded by 


E 


7 r 


2n—2 


2n(2n + l)(27r) 2n 


1 1 

7r 2 2n(2n + 1)4™ < 247r 2 

n=1 v 7 


■E4- 

71—0 


and (13.35b follows. 

Put f(Q) := 01og2 — CI 2 (tt — 9). Then f'(9) = — logsin((7r — 9)/ 2) ^ 0 and so f(9) is increasing on [0, 7r) 


and therefore f (9) ^ f( 0) = 0. provine 13.36b. 

Lemma 3.8. For 1 ^ m < k/2h zve have 

□ 

/ 7 \ 771 

(3.37) 

Ylf!(h/k)<c(h)2 k /( 2h )- m . 

(3.38) 


Proof. From Proposition l3.6l and 13.35b . 


II m(h/k) < c(h) exp (m ^1 — log(27r mh/k) + (2 tt mh/k) 2 /54^ 

< c{h) exp (m (l + log(/c/ (2nmh)) + 7r 2 /54)) 
k 


= c{h) 


— e 71 ’ 2 / 54 - 

27t mh 


< c(h ) 


te) ■ 


Similarly, Proposition 13.61 and 13.36b give 13.38b . 


□ 



Define g{x) := ilog(l/x). Then g is an increasing function on [0,1/e] with g( 0) = 0 and g( 1/e) = 1/e. 

Proposition 3.9. Let W > 0. For 5 satisfying 0 < S ^ 1/eand 61og(l/<5) < W zve have 

__ 1 , , , , N / kW\ „ mh „ 1 mh 

Il m (h/k)^ c (h)e xp l — ) for 0 < — < <5 and --$< — < 1. 

Proof. The result is true for m = 0 since J| Q l (J l /k) = 1 < c(h). For 0 < ^ ^ S, starting with 13.37b . 


n:. , (w«c W (A)' 


= c(h) exp ( ^5 


< c{h) exp ( -g (6) 


< c(h) exp ( ~W 
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since g(6) < W. For \ < ^ < 1 we have already seen in 113.341 that Y[,n(h/k) < c(h ). For 5 — <5 < ^ < 5 
we have, starting with (13.38b . 


rim^/fc) < c(h) exp 
^ c(h) exp 
^ c(h) exp 
^ c(h) exp 

Lemma 3.10. For 1 < L and 1 < m < k/2h we have 


( l0g2 2M'- 

( log2 i' 2 ’ ri / 

(r, g(S> ) 

O)- 


2nmh 

k 


□ 


\T L (m,h/k)\ < 

Proof. Combine the bounds (13.251 . (13.261 and use 

h/k 


2L-1 
2tt em 


2L — 1 




1 — ? 7 i h/k m(k/(mh) — 1 ) m 


□ 


3.4 Controlling the error term 

We have seen with Proposition 13.11 that for 1 ^ m < k/h 

n " ( '* /t) = Usm(L,/fc) ) eXP (i Cla ( 2 ™' , /‘)) 

x “P Y. (§y? (x) cot' 2 '- 2 ’ exp(-T L (m,h/k)). (3.39) 

For /. large, what is the effect of removing the factor exp (—Tj,(m, h/k)) above? Our bound on 1).(m. h/k) 
from Lemma r3.10l is poor for rn small but gets much better when m > (2 L — 1 )/ (27re). Proposition |3.9l proves 
that Yl,n(h/k) is small for small m, so we may assume m > 6k/h for a fixed 6 > 0 . As m increases we have 
n ,;, 1 (h/k) getting bigger and the bound from Lemma 13.101 getting smaller. Our goal in this subsection is 
to choose an integer L (depending on a large parameter s where 0 < h < k C s) so that these competing 
bounds produce a small enough error. 

We first give some preliminary results that shall be required in Proposition l3.12l below. Given A > 0, we 
will need real numbers r £ [ 0 , 1 ] that satisfy both the inequalities 

^( 1 + rl 0 g i) ^ Alog K’ (3 ‘ 40) 

ffP+l < 2^eA. (3.41) 

Note that the left side of (13.401 is a decreasing continuous function of r £ [0,1], decreasing from the value 
1/e to become negative at r = 1 . The left side of (13.411 is increasing and continuous from 0 at r = 0 to 1/e 2 
at r = 1. Therefore, there exist n = ri(A), r 2 = r 2 (A) so that 

j/_(l + ri logii)=Alogi (3.42) 

J&T = P-«) 


where we assume 

A < —- w 0.0079 

2 , Ke 6 
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A 

A log 1/A 

n(A) 

r 2 (A) 

Ra 

0.0079 

0.0382 

0.276 

0.924 

51.9 

0.007 

0.0347 

0.282 

0.581 

72.6 

0.006 

0.0307 

0.288 

0.427 

130.7 

0.005 

0.0265 

0.295 

0.320 

665.2 

0.00477 

0.0255 

0.297 

0.298 

11701.6 


Table 1: Some values for A and related quantities. 


so that (13.43b has a solution. If n sj r 2 then the set of all r £ [0,1] satisfying both (13.40b and (13.411 is the 
interval [n, r 2 ]. Calculations displayed in Table[l]show that ri(A) < r 2 (A) for 0.0048 ^ A ^ 0.0079. 

Suppose u > 0 and L — 1/2 = 7reA • u. If 0.0048 ^ A ^ 0.0079 then we can find r satisfying both 13.4011 
and 13.411 such that 

2i - 1 = ^1“ {3AA) 

since we may take r = r 2 (A). We next show that 13.441 still has a solution r £ [r-[. r 2 ] if we replace L by the 
integer |_7re A • uj. All this requires is that r\ and r 2 are far enough apart. Define Ra as 

r a := 3/ - ^r) (n<r 2 ). 

Lemma 3.11. Given A satisfying 0.0048 ^ A ^ 0.0079, suppose u £ R satisfies u ^ Ra- Set L := [neA ■ u\. 
Then 

2L-1 = -^u (3.45) 

for r satisfying both 13.401 and 13.411 . 

Proof. Since, as we have seen, r/e r+1 increases from 0 to 1/e 2 with r £ [0,1]/ we may use 13.451 to define r. 
From the definitions of L and r 2 we obtain 


r 2 

U 2e r2+1 


— 1 < L ^ u 


r 2 

2e r2+1 


and hence 


r 2 

V 2+1 


-3<2L-Wm 


r 2 

e r2 + 1 


- 1 . 


The right inequality in 13.461 implies that r <r 2 . Also u B Ra implies 


(3.46) 


u 


r i 

e ri+1 




u 


r 2 

e r2+1 


-3 


so that the left inequality in 13.461 implies that r\ < r. Then r £ [ri, r 2 ] implies r satisfies 13.401 and 13.411 as 
required. □ 

Proposition 3.12. Suppose A and W satisfy 0.0048 < A ^ 0.0079 and A log 1/A < IT’. For the integers h, k, s 
and m we require 

0 < h < k ^ s, Ra ^ s/h , A s/h ^ m ^ k/(2h). 

Then for L := |yreA • s/h\ we have 


Uf!(h/k)T L (m, h/k) | < (tt 3 /2 )c(h) ■ e sW /\ 
\T L (m, h/k)\ < 7r 3 /2. 


(3.47) 

(3.48) 


Proof. We write 2L — 1 = Bs/h for some B = r/e r+1 and r satisfying both 13.401 and 13.411 by Lemma [3.1 II 
With the bounds from Lemmas 13. 8l and l3.10l we have 


Um( h / k ) T L(m,h/k) < c{h) 


k \ m 7T 3 / /3s \ 
mhJ 2 \2nemh J 


13s/h 


(3.49) 
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Taking the h/s power of both sides, we see that (13.47b follows if we can establish that 


/ S \ rnh/s / /3s V 

\mhJ \2nemhJ 

or equivalently, for t = s/(mh) and 2 < t ^ 1/A, that 


fP+i/t (JL 

\27re 


< e 


< e 


w 


w 


(3.50) 


We see that has maxima on the interval (0,1 /A] at t = 1 /A and t = e r+1 . To prove (13.50b we therefore 
just need to verify it at t = 1/A and t = e r+1 . 

Since (1/A) A sC e w by the definition of A and 


P 


27reA 


< 1 


(3.51) 


by 13.41b we see that 13.50b is true for t = 1/A. Next, a short calculation shows that, for t = e 


r+1 


iog G^) ) = ^( i+riog ^)- 


Therefore 13.40b implies that 13.50b is true for t = e r+l . We have proved 13.47b . 
We also have 

7T 3 / pS V S,h < 7^ 


| T L {m,h/k)\^ — l-^—) < — ( 

2 \2iremh) 2 \27reA J 2 

using 13.51b . This proves the inequality 13.48b . □ 

As a numerical check of Proposition 13.121 take for example A = 0.006, W = 0.031, s = 500 and h = 1. 
Then L = |_ neAs/h\ = 25 and we require 3 ^ m ^ k/2 ^ 500/2. For these m and k we find 


max 

m,k 


T\,n( l l k ) T ^m,l/k) 


144.7 


max |T 2 5 (m, l/k)\ ~ 0.002 


< 8.54 x 10 7 « (tt 3 /2)c(1) • e 50(W/1 , 

< 15.5 ~ 7t 3 /2. 


(3.52) 


(The maximum of the bound on the right of 13.49b is 8.22 x 10 6 , closer to the right side of 13.52b .) Similarly, 
with the same A, W and s, take h = 3 so that L = 8 and l^ra^/c/6^ 500/6. For these m and k we find 


max 

m.k 


II ,n^/ k ) T s{m,3/k) 


0.133 


max \Tg(m, 3/fc)| ~ 0.005 


< 1.4 x 10 4 « (tt 3 /2)c(3) • e 500W/3 , 

< 15.5 « tt 3 /2. 


Since the bounds for Tjfrri. h/k) used in the proof of Proposition 13.121 (coming from Lemma 13.10b are 
independent of h/k, a short verification shows the following generalization of Proposition 13. 121 needed in 
KTSal . 

Corollary 3.13. Let W, A. s, h , k, m and L he as in Proposition ^. 12\ Suppose also that 0 < u/v < h/k. Then 

U m \h/k)T L (m,u/v) < (tt 3 /2)c(/i) • e sW/h , (3.53) 

\T L (m,u/v)\ < 7t 3 /2. (3.54) 

Proposition 3.14. For W, A, s, h, k, m and L as in Proposition ^. 12\ we have 


Yl/nWk) = 


2k sm.(ivmh/k ) 


1/2 


exp 


(2^ C1 2(“A) 

/ L—l 0 / , \ 2£—l 

I 7 Tfl \ 




(A) ~ cot< 2 '- 2 ' J + O (e"*V‘) (3.55) 


for an implied constant depending only on h. 
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Proof. With (I3.39t we see that (13.55b follows if we can prove 

n „(h/k) = n 7n(h/k) exp {T L (m, h/k )) + O (e sW ' h ) . 
For any tt > 0, say, note that the simple inequality 


e K - 1 

\e x — 1| ^ \x\ - for x € (— oo, k] 


(3.56) 


follows from the fact that ( e x — l)/x is positive and increasing. Then, using Proposition 13.12l 


Yl,n( h / k )[ ex P( T L(m,h/k))-l] ^ Y\J-(h/k)T L (m,h/k) 

< (e" 3 / 2 - l) c(h)e sW ' h . 


j7T 3 /2 _ ^ 
7T 3 /2 


□ 


Remark 3.15. The requirement A s/h m in Propositions 13.121 and 13. 14l is essential since the bound we are 

using from Lemma (3.101 


\T L (m,h/k)\ < — 


7T 3 (2L- 1 


2L —1 


worsens dramatically for 
See inequality 13.51 1 . 


2 y 27 t era 
2L — 1 ft 2neAs/h > 2nem. 


4 Expressing A\(N, a) as an integral 

4.1 First results for A\ (N, a) 

Rewrite (11.22b more in terms of N/k as 

A 1 (N,a) = Im £ 2( “ 1)?C 


%-<k^N 


k 2 


■ exp N 


iir f N 
T l "fc +3 


iir N 


xexp '-y t exp U 


N' 
2*7T(J — 
k 


and define 


9e{z) ~ ”(2^! ( 7I ’^) M lcot(2£ 2) ( 7r2: )- 


) U~N- k (m ( 4 - 1 ) 


(4.2) 


Theorem 4.1 (Sine product approximation). Fix W > 0. Let A be in the range 0.0048 < A ^ 0.0079 and set 
a := Ane. Suppose S and 5' satisfy 


^ < 6 < -, 0 < 5' < - and 6 log 1/ S, S' log 1 /S' < W. 


Then for all N ^ we have 


N 


IlN-k(m = 0(e WN ) for — G [1,1 + 5] U [3/2 — S',2) 


(4.3) 


(4.4) 


and 


riiv—fc(i/^o— 


N/k 


2Asin(7r(7V/fc — 1)) 


1 / 2 f N 

ap i^m cli{2,,N,k) 

^(lllf5^r) +0 t” J e(l + «,3/2-J') (4.5) 


with L = [a ■ AJ. The implied constants in (14.4b . d4.5b are absolute. 
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Proof. The bound d4.4b follows directly from Proposition 13.91 vvi th m = N — k and h = 1. Next, in Proposition 
13.141 we set s = N and again in = N — k and h = 1. The condition A log 1/A ^ W we need for that result 
follows from 14.3b since A<A/(1 — A) < S and A log 1/A is increasing. We also see from Table|l]in Section 
13.41 that our choice of A £ [0.0048,0.0079] ensures Ra is finite. The condition on m in Proposition 13.1 41 is 
equivalent to 


1 + 


A 

1 - A 



So 14.5b follows from Proposition l3.14l if A/( 1 - A) < 5, as we assumed. 


□ 


We will later fix some of the parameters in Theorem l4.ll 

Corollary 4.2. Let W = 0.05 and a = 0 . 0067 re ~ 0.0512. Then for all N ^ 131 we have that 14.4b . 14.5b hold when 
L = [a ■ iVJ and 

0.0061 0.01. (4.6) 

It follows from 14.1b and Theorem [44] that 


A\ ( N, a) = Im 


E 


2(—1)* 


* : £e(l+S,f-<5') 
N/k 

2Asin(7r(7V/fc — 1)) 


k 2 


exp N 


Cl 2 (2'itN/k') in f N o 
2 tt N/k + T \~k + 


1/2 


in N 


exp '-TT exp \n 


i 


N] ^g t (N/k)\ 

2’ina— + ^ N2 ^ ± 

l=i ) 


+ 0( e WN ). (4.7) 


To describe this concisely we use the notation, with z £ (1,2) to begin. 


r(z) 

Q(z) 

v(z\ N, a) 


C1 2 (2t tz) ni 


2sin(7r(z — 1)) 


1/2 


2wiaz | ge(z) 


N 


t=i 


N 2£-l 


exp(— inz/2), 

(L=[a-N\). 


(4.8) 

(4.9) 

(4.10) 


(If we need to show the dependence of v(z\ N, a) on a we may write v{z: N , a, a).) For z = z(N, k ) := N/k, 
set 

M(N, a) := E ex P (N ■ r ( z))q(z) exp(v (z; N, a)) 

k : ze(l+(5,§-5') 

and 14.7b now implies that for cr 6 Z and an absolute implied constant 

A 2 (N,a) =A 1 {N,a) + 0{e WN ). (4.11) 


We wish to replace the sum defining A 2 (N, a) with an integral. Our goal in the rest of this section is to 
prove the following, (with a = 0 . 0067 re as in Corollary 14.2b . 

Theorem 4.3. For W = 0.05 and an implied constant depending only on a, we have 


2 /*l-49 

A 2 {N,a) = — 3 ^ Im J exp(N[r(z)-Tri/z})q(z)exp(v(z-N,a))dz + 0(e WN ). 

In the next subsections we develop properties of r(z), q(z) and v(z: N, a) considered as functions of 
z = x + iy in a vertical strip in C. 


4.2 Properties of v(z; N, a) and q(z) 

Lemma 4.4. Suppose (3, 7 > 0 satisfy fiy < 1. Then given any N, d ^ 1 we have 

NAi / .,j , 1 


j=d 


for an implied constant depending only on (3, 7 and d. 
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Proof. Note that, as j varies, ) J decreases for 1 < j ^ N/(e 7 ) and increases for j ^ iV/(e 7 ). Consequently, 
for d + 1 < j < [N-/3\, 


^ max ■ 


7 (d + m (d+1) 


N 


m 


LJV./3J ■ 


<C max | N d \ (/Ly) j 
<C fV-d - 1 

since /3j < 1 implies (/Ty)^ decays exponentially with N. Thus 

d 


L^v-/9J / -w 
7 


E 


N 


<C 


7 c? 

iV 


(Ljv- jsj -d) 


1 


<c 


1 


Nd+1 ]gfd 


as we wanted. 

Proposition 4.5. Suppose 1/2 < Re(z) < 3/2 and \z — 1| ^ e > 0. AZso assume that 

f, 1 7re 

max 1 + 16 > < —. 

I e J a 

Then, for an implied constant depending only on e, a and d, 

E « wh=t (0 2,L = L«-ivJ). 




□ 


(4.12) 


(4.13) 


Proof. We have 


fff(-) <c I # 2*1 f n\z\ 


N 2 *-1 (2!)! V iV 

21—1 


2l-\ 


cot^ 2f 2 \7tz) 


<C 


J£L 

27V 


cot*- 2 ^ 2 ) ( 7 Tz) 


using J3.18I) . Suppose f ^ 2 and write z = 1 + w. Then (13.15b from Theorem 13.31 and j.3.23b show 

2C— 1 


9t{z) 

N 2(-l 


<C 


2ireN 


1 


\w 


2£-l 


*21-1 \ e -vr|y| _ 


(4.14) 


Since |z|/|w/| ^ (1 + |k/|)/|k/| < 1 + 1 je and also 8\z\ < 16 if \y\ ^ 1, we see that 


L—l 


2 L 


E 


l=d 


£=2d-l 


l + \ l 
2ne N 


£ 2 L 


V I — — 

^ \ 2TreiV 


£—2d—l 


and the proposition follows in this case with an application of Lemma [4.41 using assumption d4.12b . The 
\y\ ^ 1 case is similar, employing (13.16b from Theorem 13.31 □ 

It is now convenient to fix the choice of constants in Corollary l4.2l for the rest of this section. We note that 
condition (14. 12b in Proposition l4.5l is met for e = 0.0061 and a = 0 . 0067 re: 


1 + 16 } « max {164.9, 16} = 164.9 < 166.6=—. 

£ J a 

We have therefore shown the next result. 

Corollary 4.6. With S, S' £ [0.0061,0.01] and z £ C such that 1 + S < Re(z) < 3/2 — S' we have 

d -1 


(4.15) 


27t*ctz ^ g t {z) 

v{z\ N, a) = ^^ + 2^Jf2l=T+° 
e=i 


1 


jy2d—1 


for 2 ^ d < L = [ 0 . 0067 re • iVJ and an implied constant depending only on d. 
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Inequality (14.151 holds because the bounds obtained for gi(z) in the above proof of Proposition 14.51 are 
similar to those obtained for XL (to, h/k) in Lemma [3. 101 In particular, for z = N/k and m = N — k, (14.14b 
implies 

gt(N/k) / 2«-i \ 2 '- 1 / 16(21-1) \ 2 '- 1 

IV 2 ^ _1 \ 2tt em J \ 2neN J 

We see that condition (14.121 1 is equivalent to 


2£~l 16(2£ — 1) 

2^iem ’ 2ireN 


(4.16) 


when 2 = N/k. The requirement A s/h ^ m, see Remark |3.151 is in place in Theorem 14.11 with s = N and 
h = 1 so that AN ^ to. Hence 

2 L — 1 » 2neAN < 2item 


and 

21—1 < 2L — 1 < 16(21 — 1) < 16(2L — 1) < ^ 

2nem 2irem ’ 2neN 2ireN 

Recall that A ^ 0.0079 so that 16A < 1. Therefore AN ^ m and 16A < 1 imply (14.1611 and (14.151 . 

Proposition 4.7. The functions q(z) and v(z;N,a) are holomorphic in z for 1 < Re(z) < 3/2. In the box with 
1 + 6 ^ Re(z) ^ 3/2 — S' and —1 ^ Im( 2 ) sC 1, 


q(z), exp(v(z-,N, a)) <C 1 


for an implied constant depending only on a £ R. 
Proof. Check that for w € C, 


—tt/2 < arg(sin(7rui)) < 7 t/ 2 for 0 < Re(u>) < 1. 

Consequently, — n < arg {z/ sin(7r(2: — 1))) < 7r for 1 < Re(z) < 3/2 and so q(z) is holomorphic in this strip. 
Also v(z ; N , a) is holomorphic here since the only poles of gi{z) are at z£Z. 

Finally, q(z) is clearly bounded on the compact box, as is exp(u(z; N, a)) by Corolla ry |4. 61 □ 


4.3 Properties of r(z ) 

We defined r(z) in (14.81 for 1 < z < 2. Use (12.271 to extend it as 


" 0 ) = 


Li 2 (e 2 ™ z ) 13 tt* 


27t iz I2z ’ 

now holomorphic in the strip 1 < Re(.s) < 2. Adding a parameter j, we get 

7T lj 1 


r{z) 


2iviz . 


-Li 2 (l) + Li 2 (e 2 ™ z )-27r 2 (j + l) 


From (12.211 we obtain the identity 


(4.17) 


Li 2 (e 27r4Z ) — Li 2 (e 


z ) ., 137ri 

-7 Tl\Z — Oj — 


(1 < Re(z) < 2) 


27t*2: 27 riz 6 z 

and substituting in (14.171 produces the alternate expression, valid only for 1 < lie( 2 ) < 2, 


/ x 7 rij , 1 

r(z) 4 -- = — 7 n(z — 3) + —— 

z 2i tiz 


Li 2 (l)-Li 2 (e- 2 ^)-27r 2 (j-l) 


(4.18) 


Lemma 4.8. Consider Im(Li 2 (e 27r * z )) as a function of y e K. It is positive and decreasing for fixed x € (0,1/2) and 
negative and increasing for fixed x £ (1/2,1). 
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Proof. We have 


-f-Im(Li 2 (e 27r «)) = ImE Li 2 (e 27ri2 )) = 27rarg(l - e 2 ” 2 ). 
ay ay 


Clearly this derivative is negative for x £ (0,1/2) and positive for x £ (1/2,1). Also, we have 

lim Im(Li 2 (e 27Hz )) = Im(Li 2 (0)) = 0 

y—t oo 

implying the function decreases or increases to zero. 

Lemma 4.9. For y ^ 0 we have | Li 2 (e 27r * 2 )| ^ Li 2 (l). 

Proof. With y ^ 0 we have | e 2 ™ 2 1 ^ 1 and 


□ 


Li 2 (e 2 - 2 )| = 


°° g 27 r imz 


E 

m—1 




E m 2 Ll2 ^^' 


□ 


m =1 


Theorem 4.10. The function r(z) is holomorphic for 1 < Re(,s) < 3/2. In this strip, for j £ R, 


Re ( r(z) + ^ 


Re ( r(z) + 


z 

irij 


J 2n\z\ 




27r| 2 ; | : 


x Cl 2 (27ra) + 7 t 2 |j/| 
x C1 2 (27tx) + 7r 2 |y| 


+ 2(7 + 1) 

-2 j 


(v> 0) 

(4.19) 

(y < o). 

(4.20) 


Proof. With (14.17b . we see that r{z) is actually holomorphic for all z £ C away from the vertical branch cuts 
(—zoo, n], n £ Z. Equation (14.17b implies 


Re r(z) + 


TTZJ 


= V 


Li 2 (l) — Re (Li 2 (e 27rl2 )) + 27r 2 (j + 1)\ x Im (Li 2 (e 2 ’ rl2 )) 


27r| | 


+ 


27r|x 


For y ^ 0 we have 

Im (Li 2 (e 2 ™ 2 )) ^ Im (Li 2 (e 2w “)) = Cl 2 (27r:r) 

by Lemma [4.81 Also, using Lemma [4.91 

Li 2 (l) - Re (Li 2 (e 2 ” 2 )) + 2^ 2 (j + 1) < 2Li 2 (l) + 2 t r 2 (j + 1) 

= ^ 2 [1/3+ 2(7 + 1)] 

and (14.19b follows from 14.21b . 14.2211 and 14.23b . 

Equation 14.18b implies 


(4.21) 

(4.22) 

(4.23) 


Re ( r (z) + — ) = ny - y 


Li 2 (l) — Re (Li 2 (e 27 ” 2 )) — 27r 2 (j — 1) \ x Im (Li 2 (e 27 ” 2 )) 


27r| 2 ; | : 


27t|: 


For y 0 we have 

— xlm (Li 2 (e~ 27 ™ 2 )) < — xlm (Li 2 (e _27r “)) = — x Cl 2 (— 27 rx) = xC 1 2 ( 27 tx) 
by Lemma [4.81 Then Lemma [49] shows 

Li 2 (1) - Re (Li 2 (e- 2wi2 )) - 2n 2 (j - 1) < 2Li 2 (l) - 2 t r 2 (j - 1) 


(4.24) 


(4.25) 


and writing 


= * 2 [1/3-20-1)] 




27r| 2 : | : 


2-k\'a 


we see that 14.20b follows from 14.24b - 14.27b . 


(4.26) 

(4.27) 

□ 
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4.4 Contour integrals 

Recall that l/(2i sin(7rz)) has poles exactly at 2 ; = m£Z. Each such pole is simple with residue (—l) m /(27ri). 
By the calculus of residues, see for example Il01v74l p. 300], 


b 

£(-dm*) 

k—a 


<p(z) 

2 i sin(7rz) 


dz 


(4.28) 


for <p(z) a holomorphic function and C a positively oriented closed contour surrounding the interval [a, 6] 
and not surrounding any integers outside this interval. Next, let a, b £ Z so that 0 < a < b. With a change of 
variables in (14.28b we obtain 


b i—1 ) k 

E V-vW*) 


_i f _ del _ dz 

N J c 2 i sin(j tN/z) 


(4.29) 


for C now surrounding {N/k \ a ^ k ^ b}. 

With (14.291 , we have 

MN,a) = exp(IV • r (z)) 2is ^ N/z) exp {v{z\N,<j)) dz (4.30) 

where we may take C to be a positively oriented rectangle with left and right vertical sides 

C L :={l+S + iy : \y\ < 1/1V 2 } , C R := {3/2 - S' + iy : \y\ < 1/1V 2 } 

and with corresponding horizontal sides C + , C~ with imaginary parts 1/iV 2 and 1 /N 2 , respectively as 
shown in FigureH] Recall that we have some flexibility with 6,5' and are free to choose them in [0.0061,0.01]. 


l/N 2 - 



C+ 



I 

C L 


C R 

I . 

l/N 2 - 

1 

1 

^ 

1 * 

3/2 

1 + 5 

c~ 

3/2 

-S' 


Figure 4: The rectangle C = C + IJ GY U C U C R 


So that the path of integration in ( 14.3011 passes midway between the poles of 1/ sin(7rlV/ z), we require 

N N 

1 + 5 = --3/2 -6' = -- for a, b G Z, 6,8' G [0.0061,0.01], (4.31) 

o+l/2 a — 1/2 

The relation in (14.311 implies 

dS N (1 + S) 2 

db ~~(b+ 1/2) 2 ~ N 

so that changing b by 1 corresponds to changing 6 by « 1/JV. Similarly for a and S'. Thus, by adjusting S 
and S', we can ensure that (14.311 is true for N sufficiently large. 

Proposition 4.11. With 5, S' chosen as in (14.311 we have 

MN, cr) = —/ exp(iV • r (z)) exp(t;(z; N, a)) dz + 0(e WN ) 

jV 7 Jc+uc- 2tsm(nN/z) 

for W = 0.05 and an implied constant depending only on a. 


25 

































Proof. The proposition follows from (14.301 if we can show f CLUCR = 0(e WN ). Note that for z on the left 
vertical side we have \irN/z — 7r(6 + 1/2) | < 1/7V and similarly on the right vertical side. Since | sin(7r(6 + 
1/2)) | = 1 it follows that for large N 


2 i sm(irN/z) 


C 1 (z £ Cl U Cr). 


Proposition 14.71 implies 


q(z) exp(v(z; N,a)) < 1 (z £ Cl U Cr). 


Theorem l4.10l with j = 0 implies 


1 


27T 


37T 


Re (r (z)) < — ( xC1 2 (27 tx) + ) (z £ ClUC^) 


and we have, using Lemma [3.71 for example, 

Cl 2 (27rx) < 0.24 if 1 < x < 1.01, C1 2 (2tt*) < 0.05 if 1.49 < x ^ 1.5. 

Therefore 


1 


27T 


37T 


Re (r (z)) < — ( 1.01 x 0.24 + — ) < 0.05 (z £ C L , N > 25). 


We obtain (14.35b for z £ Cr in the same way. Consequently 

exp ( N ■ r ( 2 )) <C exp(0.051V) (z £ Cl U Cr). 

The proposition now follows from the bounds (14.321 . (14.33b and (14.361 . 

We note for future reference that 

xC1 2 (2ttx)/(2tt) < 0.04 (x £ [1,1.01] U [1.49,1.5]). 

4.5 Integrating over C + and C~ 

For the integral f c + uc - in Proposition !!. ! 11 we consider separately f c+ and f c _. Now 

1 p —tviN/z 

__ = £ (Im(z) > 0) 

j< 0, odd 



2 i s'm(nN/z) 1 — 

so that 

■f° + j< 0,odd^ G+ 

Similarly, 

1 

2 i sm(nN/z) e 2 ' 

and 



p niN/z 

X r irij N/z 


T- £ 


(Im(z) < 0) 


j> 0, odd 


/ =— V' / exp(N[r(z) + nij/z])q(z)exp(v(z-,N,a))dz. 

JC- JC- 


j> 0, odd 


The contributions to (14.391 . (14.401 when |;/| > N 2 are next shown to be negligible. 
Lemma 4.12. For N ^ 2 and to ^ 0, 

m(—N 2 — 1 — m)' 


Re r (z) + ■ 


Re ( ?' (z) + 


s? -- 


N 2 


-i(N 2 + 1 + m) \ to 


^ - 


N 2 


(z £ C+) 
(zeC~). 


(4.32) 

(4.33) 


(4.34) 

(4.35) 

(4.36) 

□ 

(4.37) 

(4.38) 

(4.39) 


(4.40) 

(4.41) 

(4.42) 
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Proof. With 2 £ C + and j = —N 2 — 1 — to. Theorem 14.101 implies 


, , ni(—N 2 — l — m)\ 1 

Re r (z) + -- I ^ 


2„|z| 2 v 10 ' 2 * 2 ”) + W 


- - 2 N 2 - 2to 
3 


(4.43) 


Since 1 < \z\ 2 < 3 and xCl2(27rx) < 3 Cl2(7r/3)/2, we see that (14.431 implies (14.411 1. The proof of (14.4211 is 
similar. □ 

With Proposition l4.7l and (14.41b it follows that 


/ exp(N [r (z) + irij / z])q(z) exp(v(z’, N,a)) dz 

. ^^2 ^AA J C + 


j<—N 2 , odd 


and this last is bounded by 


<C y' / exp(iVRe [r (z) + nij/z] ) dz <C y / exp (N(—m/N 2 )) dz 

AT2JC+ .„JC+ 

i< — N 2 ra^O 


E« 

m^O 


j — m/N _ 


1 - e~ x l N 


= 0(N). 


The same is true for j > N 2 on C and therefore the total contribution to (14.39b and (14.401 1 from terms with 

|j| > N 2 is O(N). 

Proof of Theorem 14.31 With Proposition 14 . 11 1 and the above arguments we have shown 


A 2 {N,a) = - 


N 3 / 2 


Im 


/ exp (N [r (z) + irij/z]) q(z) exp(r>( 2 ; N, a)) dz 

-N 2 ^j< o, j odd ^ c+ 


y / exp (N [r (z) + irij/z}) q{z) exp(r>( 2 ; N, a)) dz 

0<j^N 2 , j odd ( ' 


+ 0{e WN ). (4.44) 


We claim that all terms in (14.4411 are O(e 0 04Ar ) except the j = — 1 term. 

Let D + be the three lines which, when added to C + , make a rectangle with top side having imaginary 
part 1. Orient D + so that it has the same starting and ending points as C + . Since the integrand in (14.44b is 
holomorphic here we see that f c+ = f D+ . We have q(z) exp(/;(z; N, a)) <C 1 for 2 £ D + by Proposition 14.71 
On the vertical sides of D + we have 

_ / , . 7T ij\ xCh(2nx) 

Re r (z) H-- <--- < 0.04 

\ Z J Z7T 

by Theorem |4d0] and (14.37b if j < — 1. On the horizontal side of D + , with y = 1, Theorem 14.101 implies 


/ 7 rij\ 1 / 3 Cl2(7r/3) 2 


■+2 u+ 1 ) 


< 0 


if j < —1. Hence, for each integer j with —N 2 j < —1, the integral in (14.44b over C + is O(e 004N ). We will 
see later that the integral with j = — 1 cannot be bounded by O(e 0 04Ar ). 

Similarly, the integral in (14.44b over C~ is O(e 0 04Ar ), this time for all odd j with 0 < j st N 2 . Hence 

A 2 (N, a) = J + exp (N [r ( 2 ) + 7ri(-l)/z]) q{z) exp(t;(z; N, a)) dz + 0(e WN ). 

We may change the path of integration from C + to [1.01,1.49]. By (14.33b . (14.36b this introduces an error of 
size 0{e WN ). □ 
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5 Asymptotics for A\(N, a) 

5.1 The saddle-point method 

Define 


p(z) 
-4 3 (iV, a) 


(r(z) 

WJ5 lm 



, = Li 2 (l)-Li 2 (e 2 ™) 

2mz 

e~ N - p{z) q(z) ■ exp (v{z\ N, a)) dz. 


(5.1) 

(5.2) 


Thenp( 2 i) is the d = 0 case of the function Pd(z) we met earlier in ( 12.301) . We have established with (14.11b and 
Theorem [43] that, for W = 0.05, 

A 1 (N,a)=A 3 (N,a) + 0(e WN ). (5.3) 

The form of ( 15.2b allows us to find its asymptotic expansion by the saddle-point method. We state a simpler 
version of jQlv74i Theorem 7.1, p. 127] that is all we need: 

Theorem 5.1 (Saddle-point method). Let V be a finite polygonal path in C with p(z), q(z) holomorphic functions 
in a neighborhood ofV. Assume p, q and V are independent of a parameter N > 0. Suppose p'(z) has a simple zero at 
a non-corner point zo € V with Re(p(,z) — p{zo)) > 0 for z € V except at z = zq. Then there exist explicit numbers 
a 2s depending on p, q, zq and V so that we have 


l dz = f]T r(» + 1 / 2 )^ + O (j^A) 


(5.4) 


for S an arbitrary positive integer and an implied constant independent of N. 

We need to set up some notation to describe the numbers a 2 S . Write the power series for p and q near zo 
as 


p(z) =p(z 0 ) +p 0 (z- z 0 ) 2 +pi(z - z 0 ) 3 -t -, (5.5) 

q(z) = qo + qi(z - z 0 ) + qz(z - z 0 ) 2 H-• (5-6) 

(We have po ^ 0 by our assumption that p'(z) has a simple zero at z (l . For simplicity we also assume that 
qo ^ 0. This corresponds to the case (p, A) = (2,1) in |Qlv74| and the case (p, a) = (2,1) in |Woj06| .) Choose 
u) £ C giving the direction of the path V through z (l : near z o, V looks like z = zq + wt for small tel 
increasing. Note that the condition Re(p(.z) — p{zo)) > 0 implies Re(w 2 po) > 0. 

We also need the partial ordinary Bell polynomials, see I Com74l p. 136], defined as 

Bij(pi,P2,P3,...)~ .. PiPzPz--- ( 5 - 7 ) 

l^l+2^2+3^3 + -=* 

^l+^2+^3H- =j 

where the sum is over all possible t\, i 2 , £ 3 , • • • G They satisfy, for example, 

00 

(pix + p 2 a; 2 -l- Y = ’^2,Bi,j(pi,p2,...)x l (5.8) 

i=j 


and are related to the usual partial Bell polynomials by Bj j(p\. p 2 , ■ ■ ■) = j\Bij(l\pi, 2 \p 2 , ■.. )/i\. The num¬ 
bers a 2s in Theorem l5 .1 I mav be found by complicated manipulations of the series (15.511 and (15.6b . see l01v74[ 
pp. 85-86,121-127]. Wojdylo in JWoj06 Theorem 1.1] found an explicit formula for them. Adapted to the 
saddle-point method, a special case of his result is 


u> 


a2s 2( W 2po) 1 / 2 ^ q2s ~ l ^ Po 

V ' 2=0 J—0 


1/2 )& 




(5.9) 
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where we must choose the square root (w 2 po) 1 ^ 2 hi (15.9b so that Re^u^po) 1 / 2 ) > 0. Note that 

_"_ = ± 1 

(w^o) 1 / 2 (po) 1/2 

so we see that the dependence of each a 2s on the path V just involves a sign, corresponding to the direction 
of the path through the saddle-point. The first cases are 


OJ 


U> 


act — 


2(w 2 p 0 ) 1 / 2 


qo, 


a 2 = 


2(w 2 p 0 ) 1/2 \Po 


92 


3 Pigi 
2 


-P29o 


Pi 


8 pi ) 


(5.10) 


agreeing with IQlv74[ p. 127]. For V, p and zg fixed and q possibly varying in (15.4b we write a 2s (q) in what 
follows. 


5.2 A path through the saddle-point 

To apply Theorem [5d] to N :i (iV. <j) in (15.2b , we need to find the saddle-point for p(z). By Theorem 12.41 the 
unique solution to p\z) = 0 for 1/2 < Re(z) < 3/2 is given by the zg we met earlier in dl 5b . 



Figure 5: The path V = V\ U V 2 U Vg through z 0 


Now we replace the path [1.01,1.49] in d5.2b with a path through zg. We noted in Section [2.31 that zg may 
be found to arbitrary precision, and later we will set the parameter v to be 

v = Jm(z 0 )/Re(z 0 ) « 0.216279, (5.11) 

but we only require for the results below that 0.21 ^ v ^ 0.22. 

Write c := 1 + iv. Let V\ be the vertical line from 1.01 to 1.01c, V 2 the line from 1.01c to 1.49c: 


V 2 = {ct | 1.01 < t < 1.49}, 


(5.12) 


and let Vg be the vertical line from 1.49c to 1.49. So the path V := V\ LJ V 2 IJ Vg goes from 1.01 to 1.49 and 
when v is given by (15.lib it passes through zg as in Figure [5] Our goal in this subsection is to prove the 
following. 


Theorem 5.2. For the path V above, passing through the saddle point zg, we have Re(p(z) — p{zg)) > 0 for z £ V 
except at z = z 0 . 


This is exactly the requirement of Theorem 15.11 and it seems apparent from Figure [6] We prove Theorem 
I5.2l bv approximating Re(p(z)) and its derivatives by the first terms in their series expansions and reducing 
the issue to a finite computation. The path V is chosen to make this argument easier and does not use the 
line of steepest descent. 

Generalizing to Pd(z), we examine Re(pd(-)) for z on the ray z = ct for c = 1 + iv with v > 0. We also 
write 

c = pe lB (0 < p, 0 < 9 < 7t/2). 

Then, using (12.16b since \e 2mz \ < 1 when Im(z) ^ 0, 


R e[p d (ct)} = Re 


—*(Li 2 (l) + 47r 2 d)e 


—i6 i e ~ id °°^ g— 2nmvt^2nimt 


2npt 


+ 


2npt 


E 


m* 


1 (— 7r(24d + 1) sin# 1 e 27vmvt sin(2nmt — 9) 
6 p np to 2 

771 = 1 / 


2 1 


(5.13) 
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Figure 6: Graph of Re[— p(z)\ for z € V 


Similarly, employing <2.3311 . <2.341 


dt 


R e\p d (ct)\ = R e[cp' d (ct)\ = ~ | Re[p d (ct)\ + ^ 


d 2 2 / 00 N 

^■R e\pd(ct)} = Re[c 2 p^(cf)] = — — ( Re[cp^(ci)] — np ^ e _2irmvt sm(2nmt + 9) 


We may bound the tails of these series: 


(5.14) 

(5.15) 


00 g—2nmvt 


E 

m—L 




n —2Livvt 




L k ( 1 - e ~ 2nvt )' 


Collecting the first L — 1 terms in (15.131 . <5.141 and <5.151 we obtain 

^Re[p d (ct)\ = R 2 (L;t) + R*(L;t) 

(with the subscript 2 indicating the second derivative) for 

7r(24d + 1) sin# L 1 


R*(L\ t) — 

A m (t) := e~ 2 ™ vt 
B m (t) := e 


+ ^ (^A m (t) cos(27 Tmt) + B m (t) sin(27T mtfj , 


+ 


1 


6 pt 3 

m= 1 

/ 2 . / 2ttp 

Vmt 2 +Sm0 V * ' rn 2 irpt 3 

- 2 ™ vt cos6 - 1 


V t 


m?irpt 3 


and 


^(L;t)| < E 2 (L-t) ^ e _ 2m , t (^3 + ^ ■ 


We see that E 2 (L; t) is a decreasing function of L and t. We have A m (t) a positive and decreasing function 
of t. Also B m (t) is a positive and decreasing function of t when t > ^^ pm - 

Let vi = 0.21 and v 2 = 0.22. Writing pie 101 = 1 + ivi and p 2 e lS2 = 1 + iv 2 we have 

1 < pi ^ p ^ p 2 , 0 < 9i < 6 < 0 2 < 7r/2. 

For w in the interval [i>i, v 2 ], we may bound A m (t), B m (t) and E 2 (L ; f) from above and below by replacing 
v, p and 9 appropriately by v :J , pj and () :] , j = 1,2. For example 

0 < A~(t) < A m (t) < A+ (t) ( v e [ui, v 2 ]) 

with 

<(t);=e - 2 ™,^ +sinfc ^ + _J_ 

and similarly write 0 < (t) ^ B m (t ) < B^(t) and 0 < E^(L\ t ) si E 2 (L; t ) ^ (L; t ). 
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Lemma 5.3. Let c = 1 + iv with 0.21 ^ v ^ 0.22. Then ^-R e[p(ct)] > 0 for t £ [1,5/4]. 
Proof. Break up [1, 5/4] into n equal segments [xj- u xf\. Then 




min R 2 {L; t) ) - E£(L; Xj- 1 ) 

te[xj-i,xj] 


(5.16) 


Let t = x* m correspond to the minimum value of cos(27rmt) for t £ [xj-i, Xj] (so that x* m equals Xj- 1 , Xj 
or a local minimum k/2m for k odd). Similarly, let t = x** m correspond to the minimum value of sin(27r?nt) 
for t £ [xj-i,Xj]. Then 


7T S1H 0 x 

min R 2 (L; t) > ---^ (aq) cos(2nmx* jjm ) + B~(xj) sm(2Trmx** m )j 

t£[Xj— i ,Xj\ Op±Xj_-y rn —i '' ' 


(5.17) 


where we must replace A m {xf) in ( 15.17b by A+(xj_i) if cos(27r?nx* m ) < 0 and replace B m (xj) in ( 15.17b by 
B+(xj-i) if sin(27r mx** m ) < 0. 

A computation using (15.1611 and (15.17b with L = 3 and n = 2 for example shows ^j-Re[p(ct)] > 0.12. □ 

We may analyze the first derivative in a similar way. We have 

e[p d (ct)\ = R\{L\t) + R\(L; t) 

for 


Ri(L ; t ) 


7r(24d + 1) sin# 


L—l 


12 pt 2 

Cm if) ■■= e~ 2nmvt 


+ ^ (^—C m (t) cos(27r mt) + D m (t) sin(27rmt)^ , 

cos 6 


m =1 

1 sin# \ 
mt m 2 2'Kpt 2 J 


D m (t) := e 


— ‘1'ixmvt 


m 2 2Trpt 2 


and 

-2-KLvt / y 1 \ 

W (L; f )Ka(L; f ):= T3 ^( s ^ + s ). 

We see that Ei(L; t ) is a decreasing function of L and t. Also C m {t) and D m {t) are positive and decreasing 
functions of t. 

Lemma 5.4. Let c = 1 + iv with 0.21 v ^ 0.22. Then ^R e[p(ct)} > 0 for t £ [5/4,3/2]. 

Proof. Break [5/4,3/2] into n equal segments and, as in the proof of Lemma 15.31 bound J^Re \p(ct)\ from 
below on each piece. Taking n = 2 and L = 3 shows ^Re[p(cf)] > 0.03 for example. □ 

Corollary 5.5. Let c = 1 + iv with 0.21 ^ v < 0.22. There is a unique solution to -^Re[p(ct)] = 0 for t £ [1,3/2] 
that we label as t 0 . We then have Re[p(ct) — p{ct 0 )\ > 0 for t £ [1,3/2] except att = t 0 . 

Proof. Check that -^Re[p(ct)] < 0 when t = 1 and ^Re[p(ct)] > 0 when t = 5/4. By Lemma [5. 3 1 we see that 
JjRe[p(cf)] is strictly increasing for t £ [1, 5/4], It necessarily has a unique zero that we label to- By Lemma 
15.41 J^R e[p(ct)\ remains > 0 for t e [5/4,3/2] . Hence Re[p(ct) — p(cto)\ is strictly decreasing on [1, t 0 ) and 
strictly increasing on (to, 3/2] as required. □ 

Before the proof of Theorem l5.21 we need a result similar to Lemma [4.81 to let us find bounds on V\ U V:>,. 

Lemma 5.6. Consider Re(Li 2 (e 27riz )) as a function ofy ^ 0. It is positive and decreasing for fixed x with |m| < 1/6. 
It is negative and increasing for fixed x with 1/4 < |m| < 3/4. 
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Proof. We have 


Noting that 


dy 


Re(Li 2 (e 27riz )) = Im(27ri log(l - e 2niz )) = 27rlog |1 - e 2 


|1 — e 27Tiz \ < 1 
1 < |1 — e 2niz \ <2 


(M < 1/6, y > 0), 

(1/4 < |z| < 3/4, y > 0) 


(5.18) 


shows that the derivative (15.18b is negative for |x| ^5 1/6 and positive for 1/4 ^ \x\ $5 3/4. Also, we have 

lim Re(Li 2 (e 2,riz )) = Re(Li 2 (0)) = 0 

y—too 


implying the function decreases or increases to zero. 

Proposition 5.7. For 0.21 < v < 0.22 we have Re[—p(z)] < 0.0 6 for z G V\ U V 3 . 
Proof. We have x fixed as 1.01 on V\ and 1.49 on V 3 . By (14.211 we know 


Re[-p(z)] 


f(y) + g(y) 

2tt\z\ 2 


□ 


for 

f(y) ■= y (Li 2 (l) - Re(Li 2 (e 2,r “))) , g(y) = xIm(Li 2 (e 27riZ )). 

If x = 1.01 or 1.49 it follows from Lemma [4.81 that g(y) is positive and decreasing. We claim that, for these x 
values, f(y) is always positive and increasing for y > 0. 

For x = 1.01, Lemma [5.61 tells us that He(Li 2 (e 2 ^ ,| '' T+ ' , h)) is a decreasing function of y. Recalling (12.231 . 
we see it decreases from tt 2 /i 2 (0.01) < 7 t 2 /6. Therefore fig) is positive and increasing for x = 1.01. 

Next let x = 1.49. We have 


4- f(v) = Lia(l) - Re(Li 2 (e 2 " iz )) - 2th/ log |1 - e 2 ™ z | 
dy 

as in (15.181 . Lemma [5.61 implies that — Re(Li 2 (e 27r,z )) ^ 0 so that 

-7-/(2/) > 7I ' 2 / 6 ~ 2 7rylog(l + e~ 2nv ) > 0. 
dy 

Since /(0) = 0, we have shown f(y) is positive and increasing for x = 1.49. 

For z G V\, so that x = 1.01 and 0 GS y Gj Y := 1.01 x 0.22 = 0.2222, 


(/ 07 2 ) + 5(0))/(2tt1.01 2 ) » 0.0558 y G [0, Y/2] 

(. f(Y) + g(Y/2))/ (27r(1.01 2 + (Y/2) 2 ) « 0.054 y G [Y/2, Y). 


For z G V 3 , so that x = 1.49 and 0 < y < Y := 1.49 x 0.22 = 0.3278, 

Re[-p(*)] < (. f(Y) + 5 (0))/(27t1.49 2 ) « 0.0462, y G [0, Y], □ 

Proof of Theorem l5.21 Let v be given by (15.111 . Then the saddle-point z ( > lies on 'P 2 , parameterized in (15.121 . 
and when t = Re(zo) we have ct = zq. Then 

= Re[cp'(cRe(zo))] = Re[cj/( 2 o)] = 0. 

t=Re(z 0 ) 

It follows from CorolIary l5.5l that \\e\jj(z) — p(zo)\ > 0 for z G V 2 and z ^ zq. We also note that Re[—p(~o)] = 
U « 0.068076. 

For zGPiU V 3 , Proposition l5. 71 implies Re[p(z) — p(zo)\ > —0.06 + 0.068 >0. □ 


^R e\p(ct)\ 
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5.3 Asymptotic expansions 

In order to apply Theorem 15. II to (15.21) we need to understand the dependence of exp(i>( 2 ; TV, a)) on TV and 
remove this dependence from the integral. The result we prove in this subsection is the following. 

Proposition 5.8. For 1.01 < Re(z) < 1.49 and |Im(z)| < 1, say, there are holomorphic functions u a j(z) and 
Cd(z', N, a) of z so that 


d -1 

exp(v(z;N,a)) = ^ 
i=o 


u<r,j( z ) , 

Ni 


(d(z ; N, a) 


for 


Q(z] N,a) = 0 



with an implied constant depending only on a and d where 1 < d < 2L — 1 and L = [0.0067re • TVJ. 

We first establish a general result. Fix M £ Z^i. Suppose we have a function / on the positive integers 
with the following property. There exist an,... clm-i £ C and K = K(M) > 0 so that 


M -1 

im - £ 


i =1 


0-2 

N 1 


I< 

JfM 


for all A £ Zj>-|. In other words 

M— l / | \ 

a«>=£^+°(a«)- 

i=l x 7 

We next show that exp(/(TV)) must have a similar expansion to (15.201 1. 
Set A := |ai| + |a 2 | 4-h \clm-i\ and 


b 3 := 


E 

i\+2i2-\ -|-(M— l)i M -i=j 


a-, a 9 .. .a 


M—l 


Zih 2 !.. 


■ lM-V- 


(5.19) 


(5.20) 


(5.21) 


Lemma 5.9. With the above M, f, a,, K, A and bj we have 


m— 1 , 

exp(/O0) ~ E W 

l=o 


e A+K 

< - 

N m 


for any m with 1 < m < M and all N £ Z^i. 
Proof. Clearly 


/(TV) = V 

l V W Z^ jyi T jyM 


for some /m(-V) with |/m(-V)| < A'. Therefore 


M—l 


exp (/(IV)) = exp £ 


Z=1 


\ //m(JV)\ 

iw) eXp (^V^J 


We have 


( M—l 


»p|E|d = [E 


TV 1-il ii! 


E 


h M -1 




\ i—1 / \ii—0 / \*M-1=0 

Note that if we replace the a^s by their absolute values in (15.211 , (15.231 we find 

/M—l 


OO 7 

ri 

TVl' 
l=o 




1=0 \ i=l 

and in particular, ll5.24b is valid for TV = 1. 


VM 

N i 


(5.22) 


(5.23) 


(5.24) 
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With Km . (Ih23t 


=»(/<«»+ (»■>(++- 

1=0 J=m \j=0 J K v 7 


Recall the inequality (13.56b 


It follows that 


e K - 1 

|e x — 1| ^ |x|- for x £ (—oo,/c], k > 0. 


exp 


V N M J 


z K 

N M K 


( f M (N) \ _ i < f M (N) e K - 1 e« - 1 


N M ' 


(5.25) 


Hence 


°o , 

E- 


j=m 

proving the lemma. 


m 


E- 


kj=0 


TV-? 


exp 


/ /m(A) \ | 

l ATM J 




<-he' ^ 

TV m 7V M TV" 


- 1 ++ x 


□ 


If (15.19b is valid for every M £ 2+ then this is an example of an asymptotic expansion. It may be written 
formally as 

OO 

n^~Y.w 

2=1 

where the right side does not necessarily converge. Lem ma l5.9l rela tes the asymptotic expansions of exp (/(TV)) 
and /(TV). See also I01v74l p. 22] for similar exercises. 

Proof of Proposition l5.8l Recall from Corol 1 ary 14.61 tha t for z £ C such that 1.01 / W.e(z) £ 1.49 we have 

d ~ l / i \ 

(5.26) 


, . 2 niaz +1 Qt{z) _. 

v{z ] N,a) = ^- + '£^i+° 

i?=i 


TV 


N^d-l 


for 2 ^ d ^ L = [ O.OOfiTre • TVJ and an implied constant, K(d), depending only on d. For j £ Z^o put 

(2 maz + gi {z )) m1 g 2 (z ) m2 9 j {z) m i 


*<7,3 


,f0) : = 


E 


m 1 +3m 2 +5m 3 +--=j 


TO i! 


TO 2 ! 


TOj! 


(5.27) 


with u CTj o = 1- Apply Lemma l5~9l with /(TV) replaced by v(z; TV, a) and a\ = 2niaz + gi(z), a 2 = 0, • • • and 
also bj = u a j(z). Let A(a, d) be a bound for |ai| + • • ■ + |ad_i| for 2 in the stated range of the proposition. Set 


d -1 


C d(z; TV, cr) := exp(u(z; TV, cr)) - ^ ■ 




,i+) 


1=0 


TVT 


Then Q(z; TV, a) is clearly holomorphic in z and Lemma f5~9l implies |Cd(~; TV, cr)| < exp(+cr, d) + K(d))/N d 
as required. □ 

5.4 Proofs of Theorems 11.41 and 11.61 

We restate Theorem ll.6l here: 

Theorem 1.6. With b 0 = 2zoe~ mz ° and explicit b\ (a). b 2 (a),... depending on a £ Z we have 


A\ (TV, <r) = Re 


,-N 


TV 2 


b i(cr) 

TV 


frm-l(g) 


o 


/HI 


-N 


v N' 


771+2 


(5.28) 


for an implied constant depending only on a and to. 
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Proof. Recall from (12.32b that e p ' Zu ' > = 1 . Proposition 15.81 implies 


A 3 (N, a) = Im 


f e- N -^.q(z).Q(z-N,a)dz 
j'=o ^ p 


(5.29) 


where the last term in (15.29b is 

1 

<C 


jV 3 / 2 




1 ‘ jyd d" ^ ]\fd+3/2 


o -ATRe(p(z 0 )) _ 


I WO | 


— N 


jyd+3/2 


by Theorem l5.2l and Propositions ^. 7l and l5. 81 Applying Theorem 15.1 1 to each integral in the first part of (15.29b 
we obtain 


I ■ «(*) ■ WO * = 2e-»'<»> ('E r (» + l/2) aa ^ + ,"; 2 j) + O 


(5.30) 


The error term in d5.30b corresponds to an error for As(N, a) of size Ofwf x /N s+J+2 ). We choose S = d 
so that this error is less than 0{\wp\~ N /N d+3 ^ 2 ) for all j ^ 0. Therefore 


A 3 (N, a) = Im 


= Im 


4 -n-p(zq) r(a + 1/2)q 2s (g ■ u a ,j) 

jyi+3/2 e Z^ 

j=0 s=0 


jys+i/2 


V iv d+3/2 y 


w. 


-N 


2d—2 

E 


min(£,d— 1) 


= Re 


0 /—j jyt+2 

t=0 s=0 

d—2 .. £ 

^0 ]yt+2 

t =0 s=0 


E b(s + l/ 2 )a 2 s (q-u <J ,t-s) 


^ ATd+3/2 J 


E A?I+2 E r ^ S + V 2 )« 2 »(9 ■ U a ,t-s ) 


o 


I u, o I 


-AT 


V _/v d+1 


Hence, recalling ( 15.3b and with 


b t (cr) := —4z ^ T(s + l/2)a 2s (g • u a , t - s ), (5.31) 

5—0 

we obtain ( 15.28b in the statement of the theorem. 

Use the formula for a o on the left of 15.1 Ob to get 

bo(cr) = —4tr(l/2)a 0 (q • u c r , 0 ) = -Uyfn ( ^ - go) (5.32) 

\2(w 2 po) 1/2 / 

which is independent of <r. The terms po and go are defined in d5.5b . d5.6b so that, using d2.34b , 

Wo -7T ie 2mz ° 2 2 ^o 

Po =P (~o)/2 =-, q 0 = q{zo) = —■ (5.33) 

zowo w 0 

The square of the term in parentheses in d5.32b is therefore 

go = ~Zp 
4po 47re 27riz ° 

We may take ui = zp since the path V 2 is a segment of the ray from the origin through zp. A numerical check 
then gives us the correct square root: 


a o(g) = 


UJ 

2(w 2 p 0 ) 1 / 2 Q ° 


iz 0 

2 ^ne™ z ° 


(5.34) 


and the formula bo = 2zoe mz ° follows. 


□ 
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N 

m = 1 

m = 2 

TO = 3 

TO = 4 

A^na) 

200 

-33.8689 

-32.5734 

-32.4829 

-32.4681 

-32.4692 

400 

2.17937 x 10 7 

2.16780 x 10 7 

2.16710 x 10 7 

2.16712 x 10 7 

2.16712 x 10 7 

600 

1.80284 x 10 12 

1.77324 x 10 12 

1.77260 x 10 12 

1.77255 x 10 12 

1.77255 x 10 12 

800 

-3.72536 x 10 18 

-3.71475 x 10 18 

-3.71444 x 10 18 

-3.71444 x 10 18 

-3.71444 x 10 18 

1000 

-2.58000 x 10 23 

-2.54119 x 10 23 

-2.54072 x 10 23 

-2.54070 x 10 23 

-2.54070 x 10 23 


Table 2: Theorem ll.6l s approximations to A\ (N. 1). 


For example, Table|2]displays how well <15.28b in Theorem 11 ,6| a pproximates A \ (N. a) for a = 1 and some 
values of m and N. 

The expansion coefficients bt(cr) may all be written explicitly in terms of w o and zq. We give h\ (a) next. 


Proposition 5.10. Each b t (<j) is a polynomial in a of degree t. For t = 1 we have 


bi(a) = 


w 0 ((2mzo) 2 


Airizn 


e niz ° 7 T ie 3niz ° V 12 


— 2irizo + 1 ■ 


(5.35) 


Proof. Note that u a j(z) is a polynomial in a of degree j. Since a 2s (q) is linear in its argument, a 2s (ciqi(z) + 
C 2 Q 2 (z)) = cia 2 s (qi( z)) + c 2 cl 2s ( 92 ( 2 )), it follows from (15.311 that b t (a) is a polynomial in a of degree t. 

With 1=1, (15.31b implies 

bi(cr) = -4*r(l/2)a 0 (<7 ■ u<r,i) - 4*r(3/2)a 2 (g • u CT ,o) 

= -4 iy/ir ■ a 0 (q) ■ u a p{z 0 ) - 2iy/ira 2 (q). 


Since 


7 TIZ 


u a ,i{z) = 2mcjz + gi{z) = — ^12cr - - + 
we see that u a ,i{zo) = mzo(12a — 1/2 + 1/wq)/6. Also (15.10b implies 


1 — e 


27 viz 


a 2 {q) = 


Qp(g) /^72_3pigi_3p 2 + 15p| 

Po V qo 2 Po qo 2 Po 8 pi 


Taking derivatives ot q 2 (z) = iz/(\ — e 2mz ) and evaluating at z = zq shows that 


qi 


= —7 Tl ■ 


1 7T i 

2zo wo ’ 


72 

7o 


7ri 2tt~ 


1 


7T7 


37T 2 


qo 'Zzo wo qo 2 2 zo wo 8 zq 2z 0 wo 2wq 

Similarly, recalling po from d5.33b and using ( 12.33b . (12.34b and their generalizations, we have 


Pi 

Po 


1 


27T7 


P2 


1 


2ni 


2tt 2 


zq 3wq ’ po 3wo Zg 3zq wo 3w$ 


Putting this all together with ao(q) from (15.34b and simplifying completes the proof. 


□ 


Assuming Theorem ll.7l - see the summary of its proof in the next section - we may now prove our main 
result. 

Proof of Theorem ll.41 Putting h/k = 0/1 in dl.16b gives 

Cou(N) = £ (^!*)(-! y-'QoiAN). 

a— 1 ' ' 

Taking the same linear combination of dl.l8b produces 

Qhka(N)=0 for N(N + l)/2 > £ 

/ t =1 ' ' h / , 


(5.36) 


h/k^^N 
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and we partition J^jv into three parts: J^loo/ -4(TV) and the rest. The sum over this third part is 0(e WN ) by 
Theorem [L7] implying that (15.361 1 breaks into 

Cou(N)+ T,(l~ 1 M- i y~ a Qhk,( N ) + Y,(l~ 1 M- i y~ a M N ^) = °(e WN ). (5.37) 

0<h/ke&wo cr=1 '■ <T < T=1 ^ 

Use (11.17b to see that 


E ( _ J (-l) /_<, Qfck«rOO = J2(e 2mh/k - lf~ j C hkj {N). (5.38) 

cr=l ^ ' j=l 

Then Il5.38b and Theorem 1 1.61 let us write Il5.37b as 


C 0 u{N)+ Y,^ ih/k - i y~ jC ^ N ) 

0<h/ke^xoo 3 — 1 


= Re 


,-N 


TV 2 


h* 

bU + -£ + ■■■ 


TV 


N n 


O 


f I w 0 1 


— N 


V _/V m + 2 


(5.39) 


for 




We claim that b* it = 0 for 0 < t < £ — 2 and = —2^oe _,rl2o (27rjzo)^ _1 - To see this, observe that 

Ai(N, a) may be replaced in Il5.37b by 43 (TV, cr), as defined in (15.2b . since (15.3b is true. The dependence of 
. 43 (IV, a) on er comes from the exp(w(^; TV, cr)) term and we have 

E _ 1) (~ 1 ) £_<T ex P(' y ( 2; ; N > ct)) = E C - l) ( _1 E CT exp(27ri2/7V) ,T exp(u(z; TV, 0)) 

cr=l ' V<T 7 cr=l k ^ 


= [exp(27r*2:/TV) — l ] C 1 exp(t)( 2 :; TV, 1)) 


/ 2wiz (2iriz)* 


V TV 2!TV 2 
= TV 731 + 


£-1 


Ulfi(z) 


Ul,l(z) 

TV 


^.lC 2 ) v ?, 2 (-) 


TV ' TV 2 

where vj 0 (z) = {2niz) e '~ 1 , and in general, employing (15.7b . d5.8b . 


vtj(z) = ^ 1/2!, 1/3!,...) • (2 niz) e ~ 1+t ■ u^. t {z). 

t =0 


Now repeating the proof of Theorem 1 1.6 1 with replaced by v\ j yields 

t 

c u = r ( s + 1 / 2 ) a 2s(q ■ vl t _ s ) 


s=0 


(5.40) 


□ 


in the statement of the theorem, with cyo = —2zoe mz ° (2mzoY 1 = b% ( _ x as desired. 

With (15.40b we may compute the coefficients c^t explicitly. For example, similar calculations to those of 
Proposition 15 .1 01 produce vy i (z) = (2nizY(6£ + 11/2 + 1/(1 — e 2,r * z ))/12 and 


(£ + l)zo(2irizoY ZoWo(27rizoY (1 £+1 £(£ + 1 ) 

Q, l =-—-h 


gTTlZa 


e 3mz 0 y6 2irizo (2Trizo) 2 J 


(5.41) 
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6 Further results 


6.1 Proof of Theorem 11.71 

To prove Theorem ll.7l we first need a general estimate for the sine product Y[ m ( h/k ), without the restriction 
0 ^ m < k/h that was in place in Section [3] Define the set 

Z(h, k) := {(/ 3 ,7) £ Z x Z : 1 ^ |/ 3 | < fc, 1 ^ 7 < k, f3h = 7 mod fc}. 

Theorem 6.1. For all m, h,k G Z with 1 < h < k, (h, k ) = 1 and 0 < m < k we have 


■log 


n -'(h/k) 


Cl2(27rm7o/i/fc) 

27t|/3o7o| 


O 


logfc 

yfk 


( 6 . 1 ) 


where (/?o, 70) is a pair in Z(h, k) zvith |/ 3 o 7 o| minimal. The implied constant in ll 6 .lt is absolute. 


Very briefly, the proof of Theorem 16.11 involves showing, with another application of Euler-Maclaurin 
summation, that 


i loe 


IE>/*) 


S(m; h, k) 


2ir 


O 


log z k 


for 


S(rn; h, k) := ^ 

{p n )ez(h,k) 


sin(27TTO7/fc) 

\h\ 


and then relating S(m\ h, k) to Clausen's integral. 

Define D(h,k) to be the above minimal value |/?o 7 o| in the statement of Theorem 16.II For example, it is 
easy to see that 

D(h,k) = l -<=>- h = ±1 mod k 

and if D(h, k) 7 ^ 1 then 

D(h, k) = 2 -<=>• h or h = ±2 mod k 

with k necessarily odd. A simple corollary to Theorcm lb.l I says there exists an absolute constant r such that 


io s| n m (v fc )l 


C1 2 (tt/ 3) log A; 

^ -—rrrr—— + T- 


27 rD(h,k) y/k ' 

The second result we need is a general bound for Qhka(N)'- 
Proposition 6.2. For 2 ^ k < N, a € 1 and s := [N/k\ 


( 6 . 2 ) 


\Qhka(N)\ < — exp ( TV 


2 + log (1 + 3fc/4) |cr 




This proposition is proved by expressing Qhka(N) as the integral of Q(z;N,a) around a small loop 
circling h/k (recall (11.1510 and bounding the absolute value of Q(z;N,a) on this loop. A refinement of 
Proposition [621 restricting the values of k to 101 ^ k ^ N, has 


\Qhk*(N)\ < — exp ( N 


2 + log (£/2 + Ck/8) , |o- 


j}) \n N - sk (h/k) 


for £ = 1.00038 and = 1.01041. Combining (16.2b with (16.3b gives 

i - (" 2+i ° gK r w8) 


(6.3) 


(6.4) 


for k ^ 101. A straightforward calculation with (16.4b then shows that, when W > Cl 2 ( 7 r/ 3 )/( 67 r) ~ 0.0538, 
we have 

Qhka(N) < e WN /k 3 (6.5) 
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for all h/k G — J^ioo except in the cases where 

h = ±1, ±2, (k ± l)/2 mod k and N/2<k^N, 

or h = ±1 mod k and N /3 < k ^ N/ 2, 


corresponding to k large and 

D(h , k) 

small. 

Hence, the three subsets of .7 ; v — (J?~ioo U A(N)) we must 

consider separately are 






C(N) := 

|/i/fc : 

N 

~2 

< k 

-a 

0 

V/ 

h = 2 or h = k — 2j, 

V(N) := 

: 

N 

~2 

< k 

-a 

0 

V/ 

, k — 1 , k+ 1 ) 

h = 2 or h — 2 j, 

£(N) := 

|/i/fc : 

N 

"3 

< k 

II 

-e 

^ l<^ 

V/ 

or h = k — 11 


with C(N), V(N) sets of simple poles of < 5 ( 2 ;; N, a) and £(N) a set of double poles. 

To describe the asymptotics of the corresponding sums of Qhkcr{N) s, recall the dilogarithm zero w 3 = 
u’fO. —1) and its associated saddle-point zq given by ([TT5j. We also need the new saddle-points 


Z 3 := 3 + log(l — w(l, —3))/(27 ri), zi := 2 + log(l — w( 0, —2))/ {2iri) 


using the notation of Section 1231 The proof of Theorem I I .61 giving the asymptotic expansion of Ai (N, a), 
extends to cover these three new cases and, for implied constants depending only on a and m, we obtain 


Qhka{ N) = Re 

h/keC(N ) 

55 Qhkcr(N) = Re 
h/kGV(N) 

55 Qhka(N) = Re 

h/ke£(N) 


w(l,-3) N ( cl{a) c m— 1 ( cr ) 

c 0 "T r “T ‘ ‘ —T 


w. 


N 2 

-N/2 

0 I 


N 


N 2 \ 

w( 0 , — 2 )~ n 
N 2 


do(N) + 


eo 


di (cr, N ) 
N 

ei(cr) 


+ ■■■ + 


N 


+ ■■■ + 


JSfm-l 

dm —1 (n, 

Nm-l 

£m— i(rr 


O 


-N 


Hb~3)| 

jym+2 

\w 0 \~ N / 2 


jV™-i 


+ 0 


+ ^ 1 ym+2 
jym+2 


( 6 . 6 ) 

(6.7) 

( 6 . 8 ) 


where N denotes N mod 2 and the expansion coefficients may be given explicitly; the first ones are 


c* 0 = -z 3 e~ mZ3 /4, do(N) = z 0 yj2e-* iz °(e-™ ° + (— 1)*), e 0 = -3 Zl e-™'/2. 

The right sides of J6.6I ). 06.71) and 06.8D are 

O(e 0 - 0357JV ), O(e 0 0341Ar ), O(e°- 0257JV ) 

respectively and so, with 116.51 ). we have completed the summary of the proof of Theorem ll.7l See IQ'Sa ] for 
further details. 

Remark 6.3. The reason that fractions h/k in J^ioo are excluded from Theorem 1 1.71 is that the bounds we 
have for the corresponding Qhko-(N) s, due to Proposition 16.21 and its refinements, are larger than the error 
we want 0{e WN ). By improving Proposition 16.21 it should be possible to reduce =^100 to just -A-\ . See also 
IQ'Sal Remark 3.6]. 


6.2 Generalizations and conjectures 

In Table[3]we give numerical evidence for Conjecture l l .51 in the case l=\ by comparing both sides of dl.l2b 
for different values of m. These entries match those of Table |2] since c \./ = —6 t (l). Compare also flO'S15l 
Table 1]. Table|4]shows the case l = 4 of Conjecture ll. 51 
The identity <1.181 ) for a = 1 says 


55 Chki(N)=0 (IV ^2) (6.9) 

h/kEJP n 
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A 

m = 1 

m = 2 

TO = 3 

TO = 4 

Con (A) 

400 

-2.17937 x 10 7 

-2.16780 x 10 7 

-2.16710 x 10 7 

-2.16712 x 10 7 

-2.16712 x 10 7 

600 

-1.80284 x 10 12 

-1.77324 x 10 12 

-1.77260 x 10 12 

-1.77255 x 10 12 

-1.77255 x 10 12 

800 

3.72536 x 10 18 

3.71475 x 10 18 

3.71444 x 10 18 

3.71444 x 10 18 

3.71444 x 10 18 

1000 

2.58000 x 10 23 

2.54119 x 10 23 

2.54072 x 10 23 

2.54070 x 10 23 

2.54070 x 10 23 


Table 3: Coniecture ll.5f s approximations to Cou(N). 


A 

TO = 1 

TO = 2 

TO = 3 

TO = 4 

Coi4(A) 

400 

-56.2851 

-58.7844 

-58.6802 

-58.6857 

-58.6545 

600 

-1.52353 x 10 7 

-1.52212 x 10 7 

-1.52136 x 10 7 

-1.52132 x 10 7 

-1.52133 x 10 7 

800 

1.44649 x 10 12 

1.47247 x 10 12 

1.47185 x 10 12 

1.47186 x 10 12 

1.47186 x 10 12 


Table 4: Coniecture ll.5f s approximations to Cm (A). 


since C/m (A) = Qhki (A - ). With the proof of Theorem 11.21 we have shown that the largest terms in 
have h/k in J^ioo and A(N). The i = 1 case of Conjecture 11.51 indicates that all the terms with h/k in .(Coo 
are relatively small except for h/k = 0/1. So we expect 

C 0 n(A)~- E °abi(N), (6.10) 

N/2<b^N 
a=± 1 mod b 

i.e. that the asymptotic expansions of both sides of (16.10b are the same. 

Can we match up the other terms of (16.911 in the same way? A clue to the asymptotics of Cm (A) comes 
from noticing how closely it matches —1 times (16.71 . This lets us expect 

C m (A)~- E C abl {N). (6.11) 

N/2<b^N, (b, 2)=1 
a=±2 — 1 mod b 

Conjecture 6.4. For the coefficients d 0 (A), di(cr, A),... of (16.7b and an implied constant depending only on the 
positive integer to, we have 


Cm (A) = -Re ( do (N ) + ~ ' + • • • + 


—N/2 


d\ (l, A) 


dm-l(l.A)' 


A 2 


A 


A” 


O 


I wo I 


-N/2 


jym+2 


Some numerical evidence for Conjecture 16.41 is given in Table [5] The to = 1 case of Conjecture 16.41 
appeared already in [Q'S15l Conj. 6.3]. 


A 

TO, = 1 

TO = 2 

TO = 3 

TO = 4 

C 121 (A) 

1000 

1001 

1.76776 x 10 y 
2.10996 x 10 9 

1.77847 x 10 9 

2.11483 x 10 9 

1.7778 x 10 9 
2.1142 x 10 9 

1.77778 x 10 y 
2.11418 x 10 9 

1.77778 x 10 9 
2.11418 x 10 9 


Table 5: Conjecture 16.41 s approximations to Cm (A). 


Continuing the pattern from (16.10b , ( 16.lib we guess that 

Ci3i(A) + C 23 i(A)~- Y, C m {N) (6.12) 

N/2<b^N, (6,3)=1 
a=±3 1 mod b 


and indeed numerical evidence seems to support this. With more work, the techniques we have developed 
to prove Theorem ll.6l and (16.7b should give the asymptotic expansion of the right side of (16.12b . 

We finally list some further interesting directions for investigation: 


40 



























(i) Rademacher's original conjecture (11.3b has been disproved, but what is the correct conjecture? As 
discussed in ||Q'S151 Sect. 6 ], for each triple hki there seems to be some initial agreement between 
Chu{ oo) and the coefficients CW (A) for small A. 

(ii) Each Rademacher coefficient Chki{N) is a linear combination of the numbers Qhka(N) for l f a < t, 
as we have seen. For a negative, on the other hand, combinations of Qhka (A) produce the restricted 
partition function: with cr = — n in Theorem 12.1 1 we obtain Sylvester's result 


N 

Piv(«)=5^[ -Qhk(-n)(N) 

k= 1 0 ^h<k, ( h,k)=l 


(6.13) 


The inner sum in brackets is Sylvester's fcth wave [Syl82], which may be denoted as Wk(N, n). See 
also »Q'S15( Sect. 4], for example. The techniques we have developed in this paper should allow 
quantification of how well (or poorly) the first waves Wi(N, n), IT/(A, n), • • • approximate pn(ti) as 
A and possibly n tend to infinity. This ties in to work of Szekeres in ISze51l who found asymptotic 
formulas for pjv(n) when n f 0.1351V 2 by using the first wave, W\ (A', n), in the decomposition (16.13b . 
He extended this result in ISze53l . removing the restriction n AO. 135 A 2 , by using a different approach 
that incorporated the saddle-point method. 


(iii) If we replace the product 11}=i 1/(1 — q :> ) in (11.2b with a different product and examine the coefficients 
in its partial fraction decomposition as the number of factors goes to infinity, can we obtain results 
similar to Theorem 11,4P For example, following Sylvester's general theory in ]Syl82[, we may replace 
the sequence 1, 2, 3, • • • with any sequence ai, 02 , 03 . • ■ ■ of possibly repeating positive integers and 
study the partial fractions of 


N 


n T 

i=i 


— q a 


(6.14) 


as A — t 00 . The coefficient of q" in (16.14b is now the number of solutions in nonnegative integers x, to 
aixi -I-+ clnXn = n, expressed by Sylvester as a sum of waves. 
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